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SUMMARY 


Frequently,  toxicologic  processes  involve  concentrations  of 
biologically  active  substances  so  low  that  to  be  detected, 
localized,  and  quantified  they  require  attachment  of  appropriate 
radiolabels  and  subsequent  autoradiography  exposure  periods  that 
may  extend  for  many  months.  The  major  objective  of  this  research 
is  to  develop  and  deliver  a  system  to  reduce  this  long  data 
collection  period  while  maintaining  the  advantage  of  spatial 
resolution  inherent  in  autoradiography. 

After  completing  concept  validation  and  after  careful  testing  of 
three  low  light  level  camera  systems  which  were  submitted  by 
vendors  for  evaluation,  a  photon-limited  camera  and  imaging 
system  by  Hamamatsu  and  a  programmable  image  display  unit  by 
Gould  were  selected  and  ordered  for  purchase.  Software  programs 
to  distinguish  and  analyze  microstructures  and  their  properties 
were  generated. 


Aeoesalon  For 

NT1S  GRAM 
DTIC  TAB 
Unannounced 
Justification 


By - , — 

Distribution/ 


Availability  Codes 


DlBt 

1 

Avail  and/or 
Spooled 

fc'1 

Page  2 


FOREWORD 


Citations  of  commercial  organizations  and  trade  names  in  this 
report  do  not  constitute  an  official  Department  of  the  Army 
endorsement  or  approval  of  the  products  or  services  of  these 
organizations . 
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L  INTRODUCTION 


The  U.S.  Army  Medical  Research  and  Development  Command  has  a  need  to 
measure  and  accurately  localize  attogram  concentrations  of  chemical  compounds  in 
animal  and  human  tissues.  Indirect  enhancement  of  the  specific  energies  of  the 
molecules  of  these  compounds  is  a  fundamental  approach  to  this  problem.59 

Enhancement  involves  the  substitution  of  elements  in  the  compound  with  their 
radioactive  isotopes  during  synthesis.  Less  energetic  radioactive  isotopes,  such  as 
14c  (0.156  MeV)  and  3h  (0.0i.861MeV),  do  not  present  difficult  specimen  handling 
problems  and  have  half-lives  of  5,570  and  12.3  years,  respectively.  They  permit 
well-designed  biological  experiments  to  be  conducted.  The  observed  mode  of  decay 
for  both  and  is  B-emission.  Depending  upon  the  isotope  being  used,  current 
methods  require  contact  autoradiographic  exposures  of  as  much  as  6  months  to  one 
year  duration  to  record  a  usable  image.90 

Scintillation  detection  and  counting  by  this  method  will  be  a  more  effective 
alternative,  for  detecting  weak  B-emitting  radionuclides,  than  autoradiography.  In 
this  method,  the  energy  of  the  B-particles  excites  a  scintillator,  which  converts  B- 
particle  kinetic  energy  into  potentially  detectable  photons.  In  this  current  effort, 
these  photonic  events  are  temporally  and  spatially  accumulated  as  a  function  of 
position.  The  problem  is  being  approached  by  demonstrating  the  feasibility  of 
rapidly  and  microscopically  imaging  low-concentration  radiolabeled  tissues  by 
developing:  (1)  methods  for  converting  the  energy  of  emitted  B-particles  to  photons 
and  (2)  working  instrumentation  for  recording  the  "csulting  scintillations  with  high 
spatial  and  temporal  resolution.^,  23,  31,  60,  91,  97 
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H.  MATERIALS  AND  METHODS 

A.  Data  Acquisition  and  Processing  Systems 

There  are  two  major  items  of  equipment  whose  purchases  are  essentia]  for  this 
contractual  effort: 

1)  The  Gould  IP8400  Programmable  Image  Display  Unit*,  and 

2)  The  Hamamatsu  C1966-20  Video-Intensified  Microscopy  System  (VIM)**. 

While  much  of  the  thorough  evaluation  of  each  of  these  items  was  carried  out 
during  the  previous  year,  it  was  completed  in  this  year  with  decision  to  purchase.  A 
request  for  approval  was  granted,  and  the  order  for  this  equipment  was  placed 
May  23,  1986. 


B.  Development  of  Computer  Software  for  Data  Analysis 

The  program  Particlet  was  obtained  from  our  USAMR1CD  collaborator,  R.K. 
Traub.  The  primary  function  of  Particlet  is  to  enumerate  and  to  calculate  the  area 
of  microstructures  having  intensity  values  above  or  below  a  user-supplied  threshold. 
The  results  displayed  to  the  user  are: 

1)  threshold  value, 

2)  particle  count,  and 

3)  total  integrated  area  of  counted  particles. 

In  this  context,  a  "particle"  is  defined  as  a  contiguous  set  of  pixels  whose  values  all 
exceed  or  equal  some  threshold,  and  which  are  surrounded  entirely  by  pixels  whose 
values  are  all  below  that  threshold. 

Two  modifications  had  to  be  made  before  Particlet  could  be  incorporated  onto 
our  DEC  VAX  11/780  system,  worked  upon  to  develop  the  data  analysis  program,  and 
delivered  to  the  Army  in  the  format  required  on  completion  of  the  contract: 

1)  The  first  required  the  alteration  of  program  code,  so  that  image  data  were 
read  in  the  form  of  byte-packed  8-bit  data  and  then  converted  to  the 
program  format  of  64-bit  data. 

2)  The  second  was  to  increase  the  system  working  disk  space  to  fifty  thousand 
pages,  because  of  the  large  virtual  memory. 


*  Gould,  Inc.,  San  Jose,  CA 

**  Photonic  Microscopy,  Inc.,  Oakbrook,  IL 
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A  program  called  Microstructural  Analysis  Routine  (MSAR)  was  then 
developed  for  this  project  by  our  computer  scientist.  This  new  program  uses 
Particlet  as  a  basic  concept  model  and  is  written  to  analyze  data  for  which  Particlet 
had  not  been  designed.  MSAR  thus  offers  a  different  data  analysis  and  output  of  the 
extremely  limited  photon  events  (Section  III.  B.). 
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ID.  RESULTS 


A.  Data  Acquisition  and  Processing  Systems 

1.  Gould  IP8400  Programmable  Image  Display  Unit 

The  Gould  IP8400  is  a  programmable  digital  image  array  processing  and 
display  unit.  It  contains  look-up  tables  and  special  hardware  processors  to 
allow  user  manipulation  of  stored  images  and  Random  Access  Memory 
(RAM)  for  image  storage.  In  its  present  configuration  it  is  attached  to  a 
VAX-11/780,  as  its  host  computer,  as  a  peripheral  device  on  the  DEC 
Unibus. 

Image  data,  acquired  from  the  Hamamatsu  C1966  VIM  processor  by 
the  VAX-11/780,  are  transferred  to  the  IP8400  RAM  at  a  minimum  rate  of 
1.5  Mbytes/sec.  Currently  the  IP8400  has  over  524  Kbytes  of  RAM  in  two 
512  x  512  x  8-bit  memory  channels  and  in  look-up  tables  and  registers. 

Two  users  can  access  the  IP8400  concurrently,  but  memory  cannot  be 
shared.  Each  display  channel  consists  of  one  memory  channel  interfaced 
to  a  Video  Output  Controller  (VOC)  capable  of  driving  a  monochrome 
monitor  or  an  RGB  color  monitor  with  composite  video  signals. 

The  master  timing  of  the  1P8400  is  selectable  and  can  be  changed  by 
software  calls  to  internal  registers.  Each  display  channel  can  be 
manipulated  by  an  associated  joystick  and  has  an  alphanumeric  overlay 
generator  for  producing  annotation  on  the  display  monitor. 

The  IP8400  is  also  equipped  with  a  high-speed  pipeline  processor,  the 
Digital  Video  Processor  (DVP),  which  will  operate  on  a  combination  of  512 
x  512  x  8-bit  memories.  It  allows  real-time  operations  such  as  multiply, 
add,  subtract,  compare,  and  reciprocal  divide. 

2.  Hamamatsu  Cl 966-20  Video-Intensified  Microscopy  System  (MM) 

This  is  an  extremely  sensitive  system  developed  for  use  in  ultra  low-light- 
level  situations.  Its  characteristics  include  a  wide  dynamic  range  (to 
image  over  different  intensity  levels),  high  spatial  resolution,  a  high 
signal-to-noise  ratio,  and  real  time  data  acquisition  capabilities. 

After  evaluating  three  low-light-level  cameras  in  the  first  year  (Zeiss 
Venus,  Dage-MTI,  and  the  Hamamatsu  Cl 966-20)  and  final  system 
evaluations  in  this  second  year,  the  Hamamatsu  Cl 966-20  VIM  was  chosen 
for  purchase  on  the  basis  of  its  sensitivity  and  photon-limited 
characteristics. 

The  system  converts  emitted  photons  into  a  proportional  number  of 
electrons  through  the  use  of  a  two-dimensional  photon  counting  tube  that 
consists  of  three  parts.  On  one  face  is  a  photocathode,  where  the  initial 
optical  image  is  focused.  The  photons  are  converted  to  electrons  here  and 
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subsequently  accelerated  to  the  second  part  of  the  counting  tube,  a  two- 
stage  microchannel  plate  (MCP),  where  they  are  amplified  by  electron 
avalanching.  The  MCP  has  two  important  qualities: 

a)  due  to  its  design,  the  spatial  integrity  of  the  image  is 
maintained  and 

b)  because  there  is  an  adjustable  gain  setting,  image  intensifier  gain 
(avalanching)  can  be  controlled  by  the  user. 

The  gain  of  the  system  can  be  operated  in  an  analog  mode,  if  sufficient 
light  is  present,  or  in  digital  mode  for  low-light-level  imaging.  Once  the 
amplified  electrons  strike  the  phosphor  screen  that  follows  the  MCP,  a 
resultant  image  is  presented  to  an  optically  coupled  low  lag  vidicon.  This 
produces  signal  pulses  which  may  then  be  discriminated,  with  an 
acceptance  threshold,  to  eliminate  thermal  noise  that  may  occur. 

The  resolution  approaches  that  of  a  Silicon  Intensified  Target  (SIT) 
camera  (300  to  350  TV  lines  as  compared  with  400  for  the  SIT).  By 
selecting  the  type  of  imensifier  and  associated  photocathode,  a  range  of 
spectral  sensitivity  from  140  to  850  nm  is  possible. 

The  following  is  a  list  of  the  selected  components: 

a)  Specially  selected  ultra-low  noise/high  resolution  VIM  photon 
counting  camera  head  with  optional  Peltier  cooling  system.  The 
Peltier  cooling  system  was  evaluated  to  determine  its  effectiveness 
in  minimization  of  thermal  noise.  In  addition,  vendor  specifications 
indicate  that  a  dark  count  range  of  5  to  10  cps  can  be  achieved  by 
carefully  selecting  the  image  intensifier  tube. 

b)  VIM  Controller  and  control  keyboard.  To  reduce  the  number  of 
incident  light  sources  and  thereby  assure  a  light-tight  laboratory 
environment,  the  keyboard  has  been  adapted  with  a  disable  switch. 

c)  DEC  DMA  (direct  memory  access)  interface  and  software.  This 
feature  directly  couples  the  C1966-20  to  the  DEC  VAX  11/780 
computer  memory,  allowing  digitized  image  data  or  control/command 
language  to  be  transferred  in  either  direction  between  the  two 
devices.  It  offers  t lie  advantage  of  reducing  the  amount  of  user 
interaction  necessary  for  operation  c;  routine  functions. 
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Development  of  Computer  Software  for  Data  Analysis 

The  Microstructrual  Analysis  Routine  (MSAR)  was  developed  to  enhance  the 
performance  of  the  USAMRICD  Program  Particlet  (Section  11. B.),  by  providing 
additional  functionality,  in  order  to  meet  the  need  to  quantitatively  assess 
images  derived  from  radionuclide-labeled  sections  obtained  by  the  System  under 
development.  MSAR  (Appendix  A  and  Appendi:-.  B)  has  the  ability  to  perform 
the  following  five  functions: 

1.  Read  Image  Data  into  VAX  11/780  .Memory 

2.  Define  Region  of  Interest  (ROl) 

3.  Select  Threshold 

4.  Determine  Adjacency  Criteria 

5.  Calculate  Microstructure  Statistics 

1.  Read  Image  Data  into  VAX  11/780  Memory 

The  user  may  read  a  monochromatic.  512  x  512  pixel  image  from  disk  into 
VAX  memory.  The  user  is  also  requested  to  supply  the  physical  dimension 
of  the  pixel  in  microns,  if  it  is  known.  This  factor  is  important  in 
calculating  the  areas  and  perimeters  of  the  microstructures. 

2.  Define  Region  of  Interest  (RQI) 

A  region  of  interest  (ROD,  which  is  defined  by  the  user,  is  a  subset  of  the 
full  field  image  on  which  microstructure  analysis  is  performed.  If  no 
region  of  interest  is  specified,  the  entire  image  field  is  used.  The  user 
selects  an  ROI  by  entering  the  X,  Y  coordinate  values  for  the  upper  left 
and  lower  right  corners  of  the  region. 

3.  Select  Threshold 


Three  types  of  thresholding  specifications  are  available  for  user  selection: 

a)  Minimum  Threshold  —  all  pixel  values  BELOW  this  value 
will  be  ignored  during  analysis, 

b)  Maximum  Threshold  —  all  pixel  values  ABOVE  this  value 
will  be  ignored  during  analysis, 

c)  Threshold  Range  —  a  combination  of  the  above  specifications  such 
that  only  pixel  values  lying  BETWEEN  the  minimum  AND  maximum, 
inclusively,  will  be  used  in  analysis. 


These  selected  thresholds  must  be  values  between  0  and  255  to  be  valid.  In 
addition,  the  user  is  given  the  option  of  specifying  a  threshold  increment 
for  both  the  minimum  and  maximum  threshold,  either  individually  or  together. 
This  allows  analysis  to  be  performed  in  sequential  steps  automatically, 
thereby  eliminating  user  interaction  for  repetitive  operations. 
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4.  Determine  Adiacencv  Criteria 


The  user  is  given  two  choices  in  determining  the  criteria  for  contiguity  of 
adjacent  points:  ].  Orthogonal:  2.  Orthogonal  and  diagonal.  If  no 
adjacency  mode  is  selected,  the  default  is  to  number  two. 

5.  Calculate  Microstructure  Statistics 


MSAR  performs  n  variety  of  microstructure  calculations.  For  each 
specified  threshold  range  of  an  image,  the  following  are  calculated: 

a)  Total  area  of  all  microstructures  (in  pixels) 

b)  Number  of  microstructures 

c)  Average  object  area 

d)  Average  X,Y  center  of  region  of  interest 

e)  Average  X.Y  weighted  center  of  region  of  interest 

f)  Total  number  cf  perimeter  points 

g)  Average  number  of  perimeter  points  per  object 

h)  Total  perimeter  length 

i)  Average  perimeter  length 

j)  Average  horizontal  and  vertical  feret  diameters 

k)  Cumulative  and  average  (per  object)  integrated 
optical  densities 

l)  Average  maximum  optical  density  (per  object) 

The  perimeter  is  calculated  in  two  ways.  The  first,  the  total  number  of 
perimeter  points,  is  the  sum  of  all  of  the  points  on  the  boundary  of  the 
microstructure.  The  second,  the  total  perimeter  length,  is  the  distance 
travelled,  in  pixels  (or  microns  if  the  user  supplied  a  pixel  to  micron 
conversion  factor),  of  one  complete  trace  of  the  perimeter.  Once  the 
microstructures  have  been  uniquely  labeled,  the  following  statistical 
calculations  are  performed  for  each  microstructure  within  a  particular 
thresho'd  range: 

a)  X,Y  location  of  the  geometric  center 

b)  X,Y  weighted  center 

c)  Area 

d)  Total  number  of  perimeter  points  and  length 

e)  Maximum  horizontal  and  vertical  diameters 

f)  Integrated  optical  density  (pixel  values) 

g)  Average  integrated  optical  density  (pixel  values) 

h)  Maximum  optical  density  (pixel  values) 

Many  of  the  calculation  results  from  the  preceding  lists  are  given  in  units 
of  length.  For  this  reason,  the  pixel  size  is  requested  at  the  onset  of  the 
program.  The  results  obtained  will  be  more  meaningful  when  MSAR  is 
allowed  to  make  calculations  using  actual  dimensions. 
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The  hardcopy  results  obtained  upon  the  completion  of  Particlet  and  MSAR 
processing  of  the  Experimental  Image  #ED0030  (Figure  1.)  are  given  in 
Tables  1  and  2.  While  each  routine  offers  the  number  of  microstructures  and 
the  average  area  of  these  within  a  threshold  range  (or  "slice  level"),  it  is 
evident  from  this  data  that  the  results  from  MSAR  are  more  extensive  and 
useful.  The  information  concerning  averages  is  helpful  in  determining  general 
data  trends.  For  instance,  the  average  weighted  X  and  Y  center  can  direct  the 
user  to  the  region  of  highest  intensity,  which  corresponds  to  the  region  with 
the  highest  concentration  of  radiolabeled  substance. 

An  additional  level  of  detail  in  MSAR  is  illustrated  in  Table  3.  Here 
information  about  specific  microstructures  is  available  for  each  unique 
microstructure  in  the  image. 
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FIGURE  1.  Concept  Validation  Image  //ED0030  for  Signal  Dosimetry 
Experiment  with  One-second  Summation  Time. 
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TABLE  1.  Tabulated  Results  from  the  Particlet  Processing  of 
Experimental  Image  #ED0030 

EDD030.  DAT 


slice  level 
10.  00000 
20.  00000 
30.  00000 
40.  00000 
50.  00000 
60.  00000 
70.  00000 
BO.  00000 
90.  00000 
100.  0000 
110.  0000 
120.  0000 
130.  0000 
140.  0000 
150.  0000 
160.  0000 
170.  0000 
180.  0000 
190. 0000 
200. 0000 
210.  0000 
220.  0000 
230.  0000 
240.  0000 
250.  0000 


numb  e r  <  area 
1 
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16 

4 
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2 

2 

2 

1 

0 

0 

0 

0 

0 


of  particles 
24B303 
1246 
47 
15 
12 
12 
12 
1  1 
10 
10 
e 
6 
6 
5 
5 
4 
4 
4 
2 
1 
0 
0 
0 
0 
0 


TABLE  2.  Tabulated  Results  from  the  MSAR  Processing  of 
Experimental  Image  #ED0030 


FILENAME  »  ED0030  DAT 
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234  40 
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THRESHOLD 
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total 
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POINTS 
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POINTS 
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total 

PEP  1  me tcp 
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F£PET 
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10D 
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10D 

PER  OBJECT 
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Ob*4  1 

RER  OBJECT 

20-253 

207 

0  03 

3*6  00 

t . 

00  2  44 

1  «3 

31036 

84  46 

23  80 
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TABLE  3.  Tabulated  Results  for  a  Particular  Threshold  Range  from 

MSAR  Processing  of  Experimental  Image  #ED0030  (Abridged) 
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To  examine  the  integrity  and  accuracy  of  the  Microstructural  Analysis 
Routine,  a  test  image  was  synthesized  (Figure  2). 

FIGURE  2.  Microstructural  Analysis  Routine  Test  Image 
Figure  2a.  Digital  Representation 
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FIGURE  2b.  Analog  Representation  from  Monitor 


This  test  image  consists  of  a  six  by  six  matrix  containing  multiples  of  six 
different  numerical  (pixel)  values  ranging  from  zero  to  fifty  in  increments  of 
ten.  These  values  were  arranged  so  that,  as  MSAR  processed  the  image  (first 
by  analyzing  the  0  to  255  threshold  range,  then  10  to  255,  and  so  on,  until  the 
lower  threshold  value  was  equal  to  fifty),  the  pixels  that  fell  within  the 
threshold  range  formed  test  patterns  representing  the  various  conditions  of 
pixel  contiguousness  (Figure  3).  Table  4  shows  the  results  of  MSAR  execution 
on  the  test  matrix  image. 
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TABLE  4.  Tabulated  Results  from  MSAR  Processing  of 
Microstructure  Detection  Test  Image 
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C.  Literature  Review 

Emphasis  has  been  placed  on  generating  and  updating  a  cumulative 
bibliography  pertaining  to  this  project  (page  27).  For  literature  to  be  included 
in  the  bibliography,  it  must  first  satisfy  one  or  more  of  the  following  criteria: 

1.  Alternative  and  comparable  approaches  to  autoradiography,  such  as  high¬ 
speed  autoradiography  and  scintillation  autoradiography. 

2.  Physical  modelling,  such  as  the  energy  transfer  between  a  radionuclide 
and  scintillant,  and  single  photon  detection  and  capture. 

3.  Specimen  preparation  techniques,  such  as  animal  tissue  and  glass-slide 
preparation,  availability  and  assortment  of  radionuclides  and  scintillants, 
and  techniques  of  scintillation  counting. 

4.  Video  intensified  microscopy,  such  as  photodetectors  and  their  available 
technology  (photomultiplier  tubes,  charge-coupled  devices,  microchannel 
plates),  signal  processing  (noise  analysis  and  filtering),  image  processing 
and  enhancement,  and  the  analysis,  display  and  communications  of  data. 
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IV.  DISCUSSION 


fn  ,/he  r'llcr°structural  Analysis  Routine  (MSAR)  was  designed  and  implemented 
to  address  a  set  of  experimental  needs  for  this  work.  Specifications  required  that 
program  execution  time  be  minimized  without  additional  computing  hJrtw.ro 
that  two  types  of  output  files  be  creat£d7SF7ach  specified  threshold  range 

Upon  program  execution,  one  file  contains  an  array  of  the  uniquely  labeled 
microitruetures  and.  in  addition,  any  information  concerning  the  creation  of  the  file 
or  necessary  for  image  display.  Another  file  contains  the  calculated  statistics  for 
each  unique  microstructure.  A  third  file  is  generated  which  gives  the  total  and 
average  microstructure  statistics  for  the  threshold  ranges  of  a  particular  image. 
Additional  specifications  pertaining  to  the  analyses  performed  by  MSAR  were  also 
stipuifitGCi 

1.  Read  Image  Data  into  VAX  11/780  Memory 

2.  Define  Region  of  Interest  (ROI) 

3.  Select  Threshold 

4.  Determine  Adjacency  Criteria 

5.  Calculate  Microstructure  Statistics 

These  are  described  in  detail  in  Section  III.  B. 

,  .  A  tast  image  staining  various  patterns,  created  by  varying  pixel  intensity 

levels  within  a  matrix,  was  produced  to  examine  the  accuracy  of  the  Microstructu-al 
Analysis  Routine.  These  patterns  simulated  microstructures  with  edge  points  on  the 
matrix  boundary,  with  diagonal  connections,  and  with  holes.  Results  indicate  that 
microstructures  are  being  detected  and  analyzed  as  designed.  ina)c^e  that 

Using  MSAR  the  linearity  of  the  system  response  as  a  function  of  integration 
time  was  assessed  (Figure  4).  In  this  plot  the  pixel  count  above  threshold  is  s^own  as 

°f  the  lnlfegratlon  time-  The  ratio  Photoevent  area  to  full  image  field 
a  percentage  and  is  printed  adjacent  to  each  data  point.  There  isga 
linearity  of  the  relationship  for  unsaturated  images  (i.e.,  when  less  than  44  5%  of 
the  image  field  is  affected)  indicating  that,  for  a  fixed  source,  the  intepat^  time 
is  proportional  to  the  cumulative  area  of  the  detected  photoevents,  untff saturation 
t  ^e7  *  c°rrelatl0n  due  t0  overlapping  events.  Saturation  begins  to  affect 
!hi  relationship  once  the  photoevent  density  reaches  a  concentration  greater 

han  44.S%  of  the  image  field  (Figure  4).  At  this  value,  continuing  accSmulafion 
results  in  a  progressively  higher  rate  of  overlapping  points.  All  points  are  counted 

pSimTsTeaS!  lead'"S  l°  a"  underesUmalit>"  <>f  Photoevents  after  the  saturation 


AREA  ABOVE  THRESHOLD  (pixel  count) 
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FIGURE  4.  Plot  of  Image  Area  above  Threshold  versus  Signal 
Integration  Time 


Page  22 


V.  CONCLUSIONS 


During  this  second  year,  final  system  evaluations  of  The  Hamamatsu  C1966-20 
VIM  Photonic  Microscopy  System  and"  the  Gould  IP8400  Programmable  Image  Display 
Unit  were  carried  out.  Upon  receiving  official  approval,  purchase  of  the  items  was 
ordered. 

A  Microstructural  Analysis  Routine  (MSAR)  was  developed  to  analyze  the 
extremely  photon-limited  events  recorded  by  the  System  under  development.  A 
program  (Particlet)  obtained  from  the  United  States  Army  Medical  Research 
Institute  of  Chemical  Defense  (USAMR1CD)  was  used  as  a  fundamental  concept  in 
successfully  accomplishing  this  task. 

Two  peer-reviewed  scientific  manuscripts  reporting  upon  significant  aspects 
of  the  work  to  date  were  prepared  and  presented,  accepted  for  publication,  or 
published  during  this  reporting  period: 

"Detection  and  Localization  of  Quantum-Limited  Events  from  Radionuclide 
Labeled  Material  by  Computer  Enhanced  Video  Microscopy."  RK  Traub,  NJ 
Pressman,  JK  Frost,  PK  Gupta,  RL  Showers,  GW  Gill,  DL  Cook  and  JK  Frost, 
Jr.  In  the  Optical  Society  of  America’s  Topical  Meeting  on  Quantum  Limited 
Imaging  and  Image  Processing  Technical  Digest,  86(6):63-66,  1986. 

"A  Quantitative  Method  for  the  Detection  and  Localization  of  Quantum- 
Limited  Events  From  Radionuclides  in  Cells  and  Tissue  Sections  by  Computer 
Enhanced  Video  Microscopy."  NJ  Pressman,  JK  Frost,  PK  Gupta,  RL  Showers, 
GW  Gill,  DL  Cook,  JK  Frost,  Jr  and  RK  Traub.  Presented  at  the  International 
Academy  of  Cytology's  Fifth  International  Conference  on  Automation  of 
Diagnostic  Cytology  and  Histology,  Brussels,  Belgium,  May  30-June  1,  1986. 
Analyt  Quant  Cytol  Histol,  (accepted  for  publication). 
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GLOSSARY 

OF  SELECTED  KEY  TERMS,  UNITS,  AND  ABBBREVIATIONS 


AChE 


attogram 
ASA 
BID  AS 


Ci 

14c 


C1966-20 


CFIAC 

cpm 

ops 

DAGE'MTI 


-Acetylcholinesterase 

-  1 8  gram 

-  American  Standards  Association,  speed  of  photographic  film 

-  The  Johns  Hopkins  Biomedical  Image  Data  Analysis  System 

-  Curie  =  3.7  x  1010  disintegrations  per  second 

-  Radioactive  Carbon  isotope  with  atomic  mass  of  14 

-  Hamamatsu  VIM  image  processing  system 

-  Cytotechnologist  Fellow,  International  Academy  of  Cytology 

-  Counts  per  minute 

-  Counts  per  second 

-  Dage-MTI,  Inc.,  Michigan  City,  IN  46360 


DAMD17-84-C-4222 


DEC 

dps 

EM 

femtogram 

feret  diameter 

FIAC 

F1TC 

flux 

GOULD 


h 

3h 

IAC 

IOD 

JHMI 

JHU 

mCi 

MCP 


-  USAMRICD  Contract  for  "The  Development  of  a  System  for 
Enhanced  Radionuclide  Detectors"  with  JHU 

-  Digital  Equipment  Corporation,  Northboro,  MA  01532 

-  Disintegrations  per  second 

-  Electron  Microscopy 

-  10~15  gram 

-  The  maximum  diameter  of  a  structure  in  a  given  direction 

-  Fellow,  International  Academy  of  Cytology 

-  Fluorescein  Isothiocyanate 

-  Rate  of  transfer  of  particles  or  energy  per  unit  area 

-  Gould  Inc.,  DeAnza  Imaging  in  Graphics  Division,  San  Jose, 
CA  95131 

-  Planck's  constant  =  6.62363  x  10"34  joule-second 

-  Isotope  of  Hydrogen  with  atomic  mass  of  3 

-  International  Academy  of  Cytology 

-  Integrated  Optical  Density 

-  The  Johns  Hopkins  Medical  Institutions 

-  The  Johns  Hopkins  University 

-  Millicurie  =  KT^  curie 

-  MicroChannel  Plate 
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M1AC 

MSAR 

O^max 

OSA 


-  Member,  International  Academy  of  Cytology 

-  The  Microstruetural  Analysis  Routine 

-  Maximum  Optical  Density 

-  Optical  Society  of  America 


Particlet 

pVa 

PI 

pixel 

PMI 

PM1AC 

RFP 

RGB 

ROI 

slice-value 

specific 

activity 

S/N 

35s 

USAMRICD 

USAMRDC 

VAX- 11/7  80 

VENUS 

VIM 

ZEISS 


-  Particle  counting  and  sizing  program  developed  by  USAMRICD 

-  Perimeter  squared  divided  by  Area 

-  Principal  Investigator 

-  Picture  element 

-  Photonic  Microscopy,  Inc.,  Oak  Brook,  IL  60521 

-  Professional  Member,  International  Academy  of  Cytology 

-  Request  for  Proposals 

-  Red/Green/Blue 

-  Region  of  Interest 

-  Boundary  value  for  a  threshold  range 

-  Activity  C.e.,  disintegrations  per  second)  in  (micro-,  milli-) 
curies  of  radioactive  species  per  unit  mass  or  concentration 

of  non-radioactive  material  (e.g.,  56  me  C^/mM  experimental 
material) 

-  Signal  to  Noise  ratio 

-  Isotope  of  Sulfur  with  atomic  mass  of  35 

-  U.  S.  Army  Medical  Research  Institute  for  Chemical  Defense 

-  U.  S.  Army  Medical  Research  and  Development  Command 

-  DEC  32-bit  computer 

-  Zeiss  low-light-level  video  camera 

-  Video  Intensified  Microscope 

-  Carl  Zeiss,  Inc.,  Thornwood,  New  York  10594 
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APPENDICES 


Appendix  A.  MSAR  Flowcharts 


Subroutine:  NFNAHE 


Supply  image  filename 

AND  READ  IMAGE  INTO 
VAX  MEMORY. 
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PSorfPT 

Fo/L 

i  MPur 


'•c'-TUM'’/ 


F.t :  u 


Re  AD 
iMPuT 
Fil?- 

HgNPGX 


/  f'f'ORE  \ 

i.  Civ, 

Nrt? 


FForf.fi 

L^tX  vjhv:^ 
wet 
ir>  csr 


Frtofi*n-  vsEit 
fr>A  fit'GicAL 
pt  WftJJIOvJ  5 

cF  Fi*F- 


O.TTR'f  ? 


sizPij,«  j 

wv  SGT 
P-AG  -  1 


P,£^ 

U‘./y,e  if^TD 

ii  L  AlU/w 

l  ffi/*  t.TZ  i'  A 


^rt.TTf  ,nFo 
HD 

TeaMij 


-H  Kov 


TV} 


Subroutine:  DEFROI 
-Alter  region  of  interest 

(ROD  OF  IMAGE  THAT  THE 
CALCULATIONS  ARE 

performed  within. 


RgTuarJ 


SUBROUTINE:  SETTHR 


start 


-Alter  threshold  values. 
These  determine  range 

OF  PIXEL  VALUES 

considered  in  analysis. 


lM£REr*QT 

crj/ 

UFPtrL  wt  =  p 


Subroutine:  ADJCRI 


S.-TOR.-T 


-Alter  rule  governing 

CRITERIA  FOR  DETERMINING 
ADJACENCY  OF 
NEIGHBORING  PIXELS. 


rETurm 
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Subroutine:  CONTI 6 
-Determine  adjacency 

OF  PIXELS  AND 
ASSIGN  PIXELS  TO 
DISTINCT 
MICROSTRUCTURES. 


STnCT 


R  rival 


Subroutine:  STAT1S 


01  n£T 


-Calculate  statistics 
on  EACH 

MICROSTRUCTURE. 


Return 


Subroutine:  ANLZST 

CRO  STRUCT  URAL 

Data  Analysis. 
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Appendix  B.  MSAR  Program  Code 

PROGRAM  MSAR 
C 

C.  .  .  Read  in  a  REAL  *4  image  from  VAX.  measure  m i c r o s tr uc t ur e 

C.  .  .  information  <e.  g.  number,  geometric  X6Y  center,  weighted  XtY  center, 

C.  ..  area,  perimeter  length,  peria-rter  count,  mailmum  horliontel  and 

C...  vertical  feret  diameters*  Integrated  Optical  Density  and  average 

C.  .  .  Optical  Density). 

C 

INCLUDE  'EDC_MSARCQM.  FOR  ' 

BYTE  INFTXT(BO) 

DATA  10PT/1/ 

C 

C.  .  .  PRINT  RUNNING  INFO  TO  USER 

C 

OPEN  (UNIT-1,  NAME-  '  CEDC.  TASKS  3EDC_MSAR_TE  XT.  FOR  ',  TYPE-  'OLD 
CREADONLY ) 

5  READ! 1.  10, End- 15). NCHAR.  (INFTXT(J), U- 1 ,  NCHAR ) 

10  FORMAT ( 0, <NCHAR>A1 ) 

T>PE  12.  ( INFTXT(J) .  J-l,  NCHAR) 

12  FORMAT ( '  '<NCHAR>A1 ) 

GOTO  5 

15  CLOSE ( UN IT— 1 » DISPOSE-' SAVE ' ) 

TYPE  16 

16  FORMAT (  ' *TYPE  CARRIAGE  RETURN  TO  CONTINUE  :  ') 

ACCEPT  17, INFTXT(l) 

17  FORMAT (All 
C 

C 

C.  .  .  DISPLAY  INFO  TO  USER  IN  MENU  SCREEN. 

C 

CALL  WRTTMP 
CALL  WRT1NF 
C 

C.  .  .  DISPLAY  MENU  TO  USER  AT  BOTTOM  OF  SCREEN 
C 

20  CALL  WRTMEN 

C 

C.  .  .  DISPLAY  MENU  SELECTION  CURSOR.  WAIT  FOR  USER  TO  SELECT  OPTION 
C 

call  SELOPT(IOPT) 

GOTO  <  100.  200,  300.  400,  500,  600),  IOPT 
GOTO  20 
C 

100  call  nfname 

GOTO  20 

200  IF  (FNMSET  .  EQ.  0)  GOTO  <700 

CALL  DEFROI 
GOTO  20 

300  CALL  ADUCRI 

GOTO  20 

400  CALL  SETTHR 

GOTO  20 

500  I F  ( FNMSET  .  EQ.  0)  GOTO  VOO 

CALL  EXCUTE 
COTO  20 

600  CALL  CLRSCR 

CALL  CLRTXT  • 

STOP 
C 

C.  .  .  OPERATION  CANNOT  BE  PERFORMED  UNLESS  FILE  HAS  BEEN  READ  INTO  MEMORY 


I  SET  NEW  FILENAME 

! CHECK  IF  FILENAME  NOT  SPECIFIED 
IDEFINE  ROI 

IADUACENCY  CRITERIA 

! SET  THRESHOLDS 

ICHECK  IF  FILENAME  NOT  SPECIFIED 
! EXECUTE 

! EX  IT  FROM  ENTIRE  PROGRAM 
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c 

900  CALL  CLRL1NC 19. 23)  !  CLEAR  LINES  19-23  (MENU  LOCATION) 

CALL  SCROLLl 19, 23)  ! SET  SCROLL  RESJON  FOR  19-23 

TYPE  910 

910  FORMAT </'  OPERATION  CANNOT  DE  PERFORMED  UNTIL  FILE  HAS  BEEN 

C  READ  INTO  MEMORY'/) 

T=SECNDS(0.  ) 

920  CONTINUE 

3  F  ( SECNDS  (  0.  >-T  .  LT.  2.0)  GOTO  920 
CALL  SCROLL ( 1 . 24 ) 

CALL  CLRLI N ( 1 9, 23 ) 

GOTO  20 
C 

C 
C 
C 
C 

c 


END 


n  o  o  o  n 
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SUBROUTINE  URTTMP 

C 

C. . .  WRITE  TEMPLATE  TO  TERMINAL 

C 

CALL  CLRTXT 
CALL  VTCURPI4. 1 ) 

TYPE  100 
TYPE  120 
TYPE  130 
TYPE  1  A  0 
TYPE  150 
TYPE  160 
TYPE  170 
TYPE  1  BO 
TYPE  1E5 
TYPE  1=0 
TYPE  200 


1  00 

FORMAT <  ' 

♦  FILENAME:  '/  ) 

120 

FORMAT {  ' 

'  ) 

130 

FORMAT (  ' 

THRESHOLD  RANGE 

ROI ' )  * 

140 

FORMAT (  ' 

LOWER: 

UPPER: 

ROWS: 

C 

COLUMNS:  '  ) 

150 

FORMAT ( ' 

THRESHOLD  INCREMENT: 

X-LEFT: 

C 

Y-TOP :  '  ) 

160 

FORMAT  {  ' 

LOWER_INCR : 

UPPER. 

_INCR: 

X-R ICHT 

C 

V 

-BOTTOM:  7 1 

170 

FORMAT  (  ' 

'  ) 

ieo 

FORMAT  <  ' 

ADJACENCY  CRITERIA: 

'  ) 

165 

FORMAT ( ' 

PIXEL  DIMENSION: 

MICRONS') 

190 

FORMAT  (  ' 

' ) 

200 

FORMAT ( ' 

PROCRAM  STATUS:  '  1 

C 

RETURN 

END 


o  n  o  o  n 
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SUBROUTINE  WRTMEN 
C 

C.  .  .  WRITE  MENU  TO  TERMINAL 
C 

CHARACTERMO  0PTI0N<6> 

INTECER  *2  OPTLEN(i>, OPTLIN1&).  OPTCOLlfc) 
BYTE  BOLD ! 4  > .  ATTOFF ! 4  1 ,  REVERS ( 4 1 
DATA  OPTION! 1)  /'FILENAME  '/ 
data  OPT  I  ON { 2 )  /'DEFINE  ROI '/ 

DATA  OPT)ON(O)  /'ADJACENCY  '/ 

DATA  OPTION! 4)  /'THRESHOLD  '/ 

DATA  OPT  I  ON ! 5 1  /'EXECUTE  '/ 

DATA  OPT  I  ON ( 6 )  /'EXIT  '/ 

DATA  OPTLEN/6.  10.  9,  9.  7.  4/ 

DATA  0PTLIN/20,  20,  20.  22.  22.  22/ 

DATA  OPTCOL/J.  20.  39,  1. 20,  39/ 

DATA  BOLD/27,  91. 49,  109/ 

DATA  ATTOFF/27,  91.  48,  109/ 

DATA  REVERS/27,  91, 55.  109/ 

C 

TYPE  90,  (BOLP<J>.  J=t.  4) 

90  FORMAT!  '<■  '4A1  ) 

DO  100  IOPT-1,6 

CALL  VTCURP(0PTLIN<IQFT), OPTCOL! 10PT) ) 
Type  95. OPTION! IOPT) 11: 10) 

95  FORMAT !  '♦ 'A10) 

100  CONTINUE 

RETURN 
END 


.•MENU  OPTIONS 


! LENGTH  OF  MENU  OPTIONS 


•ESC  C  1  m 
•ESC  t  0  m 
!  ESC  C  7  m 


non 
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c 

c 

60 


c 

c 

90 


95 


SUBROUTINE  SELOPT(IOPT) 
CHARACTER»10  OPTIONS) 

INTEGER  *2  0PTLEN<6>.  0PTLINI6), 
BYTE  B0LD<4>,  ATT0FF(4),  REVERSE 
DATA  OPTION! 1 )  /'FILENAME  '/ 
DATA  OPT  I  ON! 2  I  /'DEFINE  ROI'/ 
DATA  OPTIONO)  /'ADJACENCY  '/ 
DATA  0PTI0NI4)  /'THRESHOLD  '/ 
DATA  OPT  ION ( 5  >  /'EXECUTE  '/ 
DATA  OPT ION l 6 )  /'EXIT  '/ 

DATA  OPTLEN/B,  10.  9,  9.  7,  4/ 

DATA  0PTLIN/20.  20,  20.  22,  22.  22/ 
DATA  OPTCOL/1 . 20. 39. 1 . 20, 39/ 
DATA  BOLD/27.  91. 49,  109/ 

DATA  ATTOFF/27,  91, 4B,  109/ 

DATA  REVERS/27,  91.  55,  109/ 

DATA  BLINK/27,91,53.109/ 

JOPT*IOPT 

TYPE  90,  (ATTOFF!  J).  J=l.  4) 

TYPE  90.  (  BOLD  <  J) .  J*=l .  4  ) 

CALL  VTCURPIOPTLIN! JODT>,  OPTCOL 
TYPE  95. OPTION! JOPT M 1 : 101 


OPTCOL (6) 

),  BLINK  14) 

•MENU  OPTIONS 


! LENGTH  OF  MENU  OPTIONS 


! ESC  C  1  m 
!  ESC  C  0  m 
( ESC  C  7  m 
! ESC  l  5  m 


IRETURN  'OLD'  OPTION  TO  NORMAL 
! HOWEVER,  BOLD  THE  OPTION 
(JOPT)  ) 


TYPE  90. (REVERS! J). J=l. 4)  ! 'REVERSE '  'NEW '  OPTION 

TYPE  90,  (BLINK(J),  J=l, 4)  !  'BLINK'  'NEW'  OPTION 

FORMAT ( '♦ '4A1 ) 

CALL  VTCURP!0PTL1N( IOPT),  OPTCOL(IOPT) > 

TYPE  95, OPTION! IOPT) ( 1:  10) 

FORMAT! '♦'A10) 

CALL  VTCURP (23,  1 ) 

WAIT  FOR  USER  TO  MOVE  MENU  CURSOR 


IKEY=0 

TLAST  =  SECNDS(0.  )  ISECONDS  SINCE  LAST  INPUT 

100  I  =  I  ESC ( I  KEY ) 

IF! IKEY  .  EQ.  0)  THEN 

IF  ( SECNDS  (  0.  l-TLAST  .  CT.  3600.  )  THEN 
I0PT=6 
GOTO  999 

END  IF 
GOTO  100 

END  IF 

IF! IKEY  .  EQ.  13)  GOTO  999  !<CR>  USER  ACCEPTS  IOPT 

IF!  IKEY  .  LT.  65  .OR.  IKEY  .  CT.  6B )  GOTO  100 
JOPT* IOPT 

IF! IKEY  .  EQ.  65)  THEN  !UP  ARROW 

10PT=I0PT-3 

IF!  IOPT  .LT.  1)  I0PT-I0PT*6 

END  IF 

IF! IKEY  .  EQ.  66)  THEN  'DOWN  ARROW 

I0PT=I0PT*3 

IF! IOPT  .  CT.  6)  IOPT-IOPT-6 

END  IF 

IF! IKEY  .  EO  67)  THEN  IRICHT  ARROW 

I0PT=10PT+1 

IF!  IOPT  .  CT.  6)  I  OPT*  I 


END  IF 


o  n  o  n  o 
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I F  (  I  KEY  .  EQ.  68)  THEN  ‘LEFT  ARROW 

I  OPT  = i OPT-1 

1F(  10PT  .  LT.  1  )  10PT-6 

END  IF 
GOTO  60 
C 

999  TYPE  90,  (ATTOFFI J).  J-l- 4)  (RETURN  OPTION  TO  NORMAL 

RETURN 
END 
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c 

c. .  . 
c 


c 

c. .  . 
c 


c 

c. .  . 
c 

eo 

90 

100 

120 


c 

c. .  . 
c 

300 

320 


330 

340 


C 

c. .  . 
c 


c 

c. .  . 
c 

350 


SUBROUTINE  NFNAME 
INCLUDE  'EDC_MSARCOM.  FOR ' 

SETUP  SOME  PARAMETERS 

IRROW-512 

IRC0L=S12 

JRXTOP-1 

IRYTOP-1 

IRXBDT=512 

IRYBOT-512 

CLEAR  MENU  SECTION  OF  VT  SCREEN 

CALL  CLRLINI 19, 231  ! CLEAR  LINES  19-23  (MENU  LOCATION) 

CALL  SCROLL  <19, 23 )  !SET  SCROLL  REC10N  FOR  19-23 

PROMPT  USER  FOR  FILENAME 

DO  90  J=l.  14 
FNAMEI J)=0 
TYPE  100 

FORMAT < 'tENTER  FILENAME  AND  EXTENSION  :  ') 

ACCEPT  120.NCHAR,  < FNAME ( U) , J-l ,  NCHAR ) 

FORMAT ( Q.  13A1 I 
IF  (NCHAR  .  EQ.  0)  THEN 

CALL  SCROLL (1,24) 

CALL  CLRLIN( 19.  23) 

GOTO  999 

END  IF 

GET  PIXEL  DIMENSIONS  FROM  USER 
TYPE  320 

FORMAT ( '  ENTER  PHYSICAL  DIMENSION  OF  PIXEL  IN  MICRONS 
C  <if  0.  25V) 

TYPE  330 

FORMAT! 'tTYPE  <CR>  IF  UNKNOWN  :  ') 

ACCEPT  340. NCHAR, SIZPIX 
FORMATIQ.  F15.  5) 

IPIXFG=0  ! SET  PIXEL  FLAG-0 

IF (NCHAR  EO  0)  THEN 

IP  I XFG= 1  ! 5ET  PIXEL  FLAG  SINCE  USER  DIDN'T  KNOW  SIZE 

SIZPIX-1.  0 

END  IF 

OPEN  INPUT  FILE 

0PEN1UNIT-3. NAME-FNAME. TYPE- 'OLD READONLY,  FORM- 'UNFORMATTED ' , 
CERR-900 ) 

READ  IMAGE  INTO  MEMORY 
TYPE  350 

FORMAT (////'  READING  IMAGE  FILE  INTO  MEMORY'//) 

READ  (3)  TEMP 

CLOSE (UNIT-3,  DISPOSE- 'SAVE ' ) 

NCOLS-512 
NROWS-5 1 2 
DO  400  J-l, 512 


.'TRANSFER  TO  1*2  ARRAY 
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4  00 

c 

c. .  . 
c 


c 

999 

C 

C 

900 

910 


C 

C 

C 

C 

c 


DO  400  K-l,  512 

1MACE1K. J)-I INTI  TEMP (K,  J> ) 

WRITE  NEW  INFORMATION  TO  TERMINAL 

CALL  SCROLL ( 1 . 24 ) 

CALL  WR T INF  IWRITE  INFO  TO  SCREEN 

CALL  CLRLIN1 19,  23) 

FNMSET-1  .'SET  FLAO  TO  INDICATE  FILE  READ  SUCCESFULLY 

RETURN 


TYPE  910 

FORMAT! ’  ERROR  IN  OPENING  FILE.  PLEASE  REENTER') 

GOTO  BO 

END 
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SUBROUTINE  DEFROJ 
INCLUDE  'EDC_MSARCOM.  FOR  ' 

C 

C.  .  .  CLEAR  MENU  SECTION  OF  VT  SCREEN 
C 

CALL  CLRLIN) IS.  23)  ! CLEAR  LINES  IB-23  (MENU  LOCATION) 

CALL  SCROLL) IB, 23)  ! SET  SCROLL  RECION  FOR  IB-23 

C 

C. , .  KEYPAD  ENTRY 

C 

200  TYPE  220 

220  FORMAT ( / '*ENTER  X.Y  COORDINATES  DF  TOP  LEFT  CORNER  OF  ROI  :  ') 

ACCEPT  240.  NCHAR,  1 1.  12 
240  FORMAT ( 0, 21 4 ) 

1 F ( NCHAR  .  EO.  0)  COTO  99? 

IFIU  .  LT.  0  .OR.  12  .LT.  0)  CDTD  300 
I F  <  I  I  CT.  NCOl_S  .OR.  12  .  CT.  NROWS)  COTO  300 
C 

250  TYPE  260 

260  FORMAT)  '*ENTER  X,Y  C0RRD1NATES  OF  BOTTOM  RIGHT  CORNER  OF  ROI  :  ') 

ACCEPT  240, NCHAR, 13, 14 
IF) NCHAR  .  EO.  0)  COTO  999 

IF  (  13  .  LT.  II  .OR.  14  .  LT.  12)  GOTO  320 

IF)  13  .  GT  NCOLS  .OR.  14  .  GT.  NROWS)  GOTO  300 
C 

1RXT0P=I 1 
IRYT0P=12 
IRXB0T=I3 
IRYB0T=I4 

IRR0W=1RYB0T-1RYT0P+1 
1RC0L=1RXB0T-1RXT0P+1 
GOTO  999 
C 

C.  .  .  ERRORS  IN  ENTERING  ROI  INFORMATION 
C 

300  TYPE  310 

310  FORMAT)'  BOUNDARY  DEFINITION  ARE  OUTSIDE  OF  LEGAL  RANGE'/) 

GOTO  200 
320  TYPE  330 

330  FORMAT)'  BOTTOM  R1CHT  MUST  BE  BELOW  AND  TO  THE  RIGHT  OF  TOP  LEFT') 

COTO  250 
C 

999  CALL  SCROLL <1. 24) 

CALL  WRTINF  ! WRITE  INFO  TO  SCREEN 

CALL  CLRLIN) 18. 23) 

C 

RETURN 

END 

C 

C 

c 

c 

c 


Page  52 


SUBROUTINE  SETTHR 
INCLUDE  'EDC_MSARCOM,  FOR  ' 

BYTE  ACEPT 
C 

C.  .  .  CLEAR  MENU  SECTION  OF  VT  SCREEN 
C 

CALL  CLRLIN! 19,  23)  ! CLEAR  LINES  19-23  'MENU  LOCATION) 

CALL  SCROLL! 19, 23)  !SET  SCROLL  REGION  FOR  19-23 

C 

C.  .  .  CET  THRESHOLD  VALUES  FROM  USER 
C 

90  TYPE  100 

100  FORMAT!  'FENTER  INITIAL  LOWER, UPPER  THRESHOLD  VALUES  :  ') 

ACCEPT  1  20,  NCHAR.  1 1,  12 
120  FORMAT ! 0,  214) 

I F ( NCHAR  .  EQ.  0)  COTO  999 
I F !  I  1  .  LT.  0  .OR.  12  .  LT.  11)  GOTO  BOO 
IF!  12  .  CT.  255  .OR.  II  .  GT.  12)  CDTO  BOO 
C 

type  160 

160  FORMAT! 'FENABLE  THRESHOLD  INCREMENTATION  ?  <N/Y>  :  ') 

ACCEPT  170. NCHAR,  ACEPT 
170  F ORMAT  (  Q,  A 1  ) 

IFIACEPT  .  EQ.  *Y'>  COTO  190 


C 

C.  .  .  USER  DID  NOT  WANT  THRESHOLD  INCREMENTATION.  SET  VALUES  AND  RETURN 
C 

ITLOW=I I 
I  TUP R* 1 2 
IT1NCR=0 
ITLINC-0 
1 TUINC-0 
GOTO  999 
C 

C.  .  .  CET  INCREMENTATION  VALUES  FROM  USER 
C 

190  TYPE  200 

200  FORMAT! '*ENTER  LOWER. UPPER  THRESHOLD  INCREMENT  VALUES  :  ') 

ACCEPT  120.  NCHAR,  13.  14 
IFINCHAR  . EQ.  0)  GOTO  999 

IF!  13  .  CT.  255  OR.  14  .  CT.  255)  COTO  B20 
IF!  13  .  EQ.  0  .AND.  14  .  EQ.  0)  GOTO  840 
C 

C. . .  SET  VALUES  AND  RETURN 

C 

ITLOU-I 1 
1 TUPR  =  I  2 
ITINCR'l 
ITHNC«I3 
I  TUI  NC*=  I  4 
GOTO  999 
C 

800  TYPE  810 

810  FORMAT!/'  INVALID  THRESHOLD  VALUES.  PLEASE  RE-ENTER'/) 

GOTO  90 
C 

820  TYPE  830 

830  FORMAT!/'  THRESHOLD  INCREMENT  VALUE  CANNOT  EXCEED  255.  PLEASE 

C  RE-ENTER'/) 

GOTO  190 
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c 

040 

850 

C 

999 


C 


C 

c 

c 

c 

c 


FORMAT^/ '  INVALID  INCREMENT  VALUES.  PLEASE  RE-ENTER'/) 
COTO  190 


CALL  SCROLL  < 1  *  24  > 
CALL  WRTINF 
CALL  CLRL1NI IB.  23) 


'WRITE  INFO  TO  SCREEN 


RETURN 

END 


non 
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SUBROUTINE  ADJCRI 
INCLUDE  'EDC_MSARCOM.  FOR ' 

clear  menu  section  of  vt  screen 

CALL  CLRLINI 19.  23)  ! CLEAR  LINES  19-23  (MENU  LOCATION) 

CALL  SCROLL< 19.  23)  ! SET  SCROLL  RECION  FOR  19-23 

C 

C.  .  .  CET  USER  INPUT  AND  CHANGE  ADJACENCY  CRITERIA 
C 

90  TYPE  100 

100  FORMAT ( /  '  SELECT  CRITERIA  BY  WHICH  POINTS  ARE  CONSIDERED 

C  CONTIGUOUS'/’  1 -ORTHOGONAL  POINTS  ONLY’/’  2 -ORTHOGONAL  AND 
C  DIAGONAL  POINTS'  1 
TYPE  120 

120  FORMAT< -*£NT£R  SELECTION  :  ') 

ACCEPT  MO.  NCHAR.  1SEL 
MO  FORMAT  (0,  12) 

1 F  ( NCHAR  .  EO.  0)  GOTO  999 
1FIISEL  .  LT.  1  .OR.  ISEL  .  GT.  2)  GOTO  90 
I F ( ISEL  .  EO.  1 )  1ADJCt«0 
IF( ISEL  .  EG.  2)  IADJCT=1 
C 

999  CALL  SCROLL ( 1 . 24 ) 

CALL  URTINF  !WRITE  INFD  TO  SCREEN 

CALL  CLRL1N1 19,  23) 

C 

RETURN 

END 

C 

C 

C 

c 

C 
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SUBROUTINE  EXCUTE 
INCLUDE  '  EDC_MSARCOM.  FOR* 

C 

C.  .  .  SAVE  INITIAL  THRESHOLD  VALUES 
C 

JTLOW-ITLOW 

JTUPR-ITUPR 

NUMAVG-0 

C 

C.  .  .  CHANCE  PROCRAM  STATUS  TO  'EXECUTING 
C 

ST  ATUS=2 
CALL  WRTINF 
C 

C.  .  .  OPEN  OUTPUT  FILE  TO  CONTAIN  'AVERACE'  INFORMATION  OF  MICROSTRUCTURES 
C 

OPEN ( UNI T -2. NAME- 'MICROS.  DAT ' . TYPE- 'UNKNOWN ACCESS- 'APPEND  '  ) 

DO  10  U=l.  14 

I F  (  FNAME  ( J )  .  EQ.  '.'1  COTQ  15 
10  CONTINUE 

NCHAR-14 
GOTO  20 

lb  NCHAR-U--3 

20  UR1TE<2,  25).  (FNAME(J).  J=l,  NCHAR),  IRXTOP.  1RYT0P.  1RXB0T.  IRYBOT 

25  FORMAT!//'  FILENAME  =  '<NCHAR>A1 // '  REGION  OF  INTEREST 

CTOP  LEFT  ( X= ' 13'. Y- '13 ' ).  BOTTOM  RIGHT  ( X- ' 1 3 ' , Y» ' 1 3 ' 1  ' // > 

WRITEI2.30)  ! PUT  OUT  HEADER  OF  FIRST  SECTION 

30  FORMAT ( ' 

C  AVERACE  AVERACE'/ 

C '  NUMBER  AVERACE 

C  AVERAGE  AVERACE  WEIGHTED  WEICHTED ' / 

C’  THRESHOLD  TOTAL  OBJECTS  OBJECT 

C  X-  Y-  X-  Y-'/ 

C'  RANGE  AREA  (MICROSTRUCTURES)  AREA 

C  CENTER  CENTER  CENTER  CENTER'/) 

C 

C 

C.  .  .  CLEAR  MENU  SECTION  OF  VT  SCREEN 
C 

CALL  CLRLIN!  19.  23)  ! CLEAR  LINES  19-23  (MENU  LOCATION) 

CALL  SCROLL( 19.  23)  ! SET  SCROLL  REGION  FOR  19-23 

C 

C.  .  .  CO  THROUGH  THE  INPUT  IMAGE  FILE  IN  ARRAY  ' I MAGE ' .  THRESHOLD  USING 
C...  THE  USER-SUPPLIED  VALUES.  AND  PLACE  THE  BINARY  RESULT  IN  'THRESH'. 

C.  .  .  ONLY  BE  CONCERNED  WITH  THE  REGION  OF  INTEREST. 

C 

50  DO  100  I-IRYTOP, IRYBOT 

TYPE  99.  1 

99  FORMAT ( '  THRESHING  LOC  '13) 

DO  100  J-IRXTOP,  IRXBOT 

IF  <  IMAGE  ( J.  I  )  .  CE.  ITLOW  .AND.  IMAGEfJ.  I)  .  LE.  ITUPR)  THEN 
THRESH  (J.  I  )  - 1 

C.  .  .  «**  NOTE  -  NEXT  LINE  OF  CODE  IS  FOR  VERTICAL  FILL  OF  IMACE. 

IF  (  IMAGE  ( J.  T  — 1  )  .  LT.  ITLOW  .  OR.  IMAGE(J.I-l) 

C.  CT.  ITUPR)  THRESH!  J.  1  —  1  >  — 1 
ELSE 

THRESH!  J,  I  )  -0 

END  IF 

100  CONTINUE 
C 
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C.  .  .  GO  BACK  INTO  IMAGE  AND  SEARCH  FOR  -1  FLAG.  IF  FOUND.  PUT  THE  VALUE 

C.  .  .  OF  THE  PIXEL  DIRECTLY  VERTICAL  TO  THE  PRESENT  PIXEL  INTO  PRESENT 

C.  .  .  LOCATION. 

C 

DO  2100  I-IRYTOP. 1RYB0T 
TYPE  99,  I 

DO  2100  U=IRXTOP.  IRXBOT 
IF ( THRESH! J,  I )  ,  EQ.  -1)  THEN 
THRESH!  J,  I  )«1 
IMAGE  ( J*  1)«1MAGE(J.  1*1) 

END  IF 

2100  CONTINUE 

C 

C. . .  CALL  SUBROUTINE  'C0NT1G'  TO  CALCULATE  THE  CONTIGUOUS  STRUCTURES 

C. . .  IN  THE  USER-DEFINED  ROI  AND  USING  THE  USER-SPECIFIED  ADJACENCY 

C.  .  .  CRITERIA.  RESULT  OF  THIS  SUBROUTINE  WILL  BE  UNIQUELY  LABEL 
C.  .  .  MICROSTRUCTURES  IN  THE  ARRAY  'THRESH'. 

C 

CALL  CONTIG 
TYPE  101.  UNIQUE 

101  FORMAT!'  UNIQUE  NUMBER  -  '17) 

T  =  SECNDS  ( 0.  > 

102  CONTINUE 

1FISECNDSI0.  )-T  .  LT.  2.0)  GOTO  102 
C 

C.  .  .  CALL  SUBROUTINE  ‘STATIS*  TO  CALCULATE  THE  STATISTICS  ON  EACH 
C.  .  .  OF  THE  UN10UE  MICROSTRUCTURES  IN  THE  ARRAY  'THRESH'.  WRITE 

C.  .  .  OUTPUT  TO  OUTPUT  FILE 

C 

CALL  ST  AT  1 S ( NUMAVG  > 

C 

C.  .  .  WAS  THRESHOLDING  INCREMENTATION  ENABLED? 

C 

IF ( IT  I  NCR  .  EQ.  0)  GOTO  900 
C 

C. . .  INCREMENT  THRESHOLDS 

C 

ITLOW=ITLOU+ITLINC 

ITUPRUTUPRUTUINC 

IFUTLOW  .  LT.  0  .OR.  ITLOW  .  GT.  255)  COTO  900  IOUTSIDE  VALID  RANCE 
IF (  I TUPR  .  LT.  0  OR.  1 TUPR  .  GT.  255)  GOTO  900 
IFUTLOW  .  GT.  ITUPR)  GOTO  900 
C 

C.  .  .  WRITE  INFORMATION  TO  TERMINAL 

C 

CALL  SCROLL ( 1 .  24) 

CALL  WRTINF 
CALL  CLRLIN! 19.  20) 

CALL  SCROLL! 19,  23) 

COTO  50 
C 

C. . .  RESET  INITIAL  THRESHOLD  AND  THEN  RETURN 

C 

900  ITLOW«JTLOW 

ITUPR«JTUPR 
C 

C.  .  .  PUT  A  <FF>  IN  FILE  TO  SEPERATE  THIS  RUN 
C 

WR1TEI2. 910), 12 
910  FORMAT!'  '.Al) 


'.CLEAR  LINES  19-23  (MENU  LOCATION) 
! SET  SCROLL  REGION  FOR  19-23 


u  o  u  u 
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C.  .  .  WRITE  SECOND  PACE  OF  INFO 
WR  I TE  (  2.  915) 

915  FORMAT ( '  AVER ACE 

C  AVERAGE  AVERAGE'/ 

C'  TOTAL  PERIMETER  A 

CVERAGE  HORIZ  VERTICAL  AVERAGE 

C  AVERAGE'/ 

C'  THRESHOLD  PERIMETER  POINTS  TOTAL  PERI 

C METER  FERET  FERET  CUMULATIVE  IOD 

C  ODma i  '  / 

C'  RANCE  POINTS  PER  OBJECT  PERIMETER  PER  0 

OBJECT  DIAMETER  DIAMETER  IOD  PER  OBJECT 

C  PER  OBU.'CT'/) 

DO  930  U-l.NUMAVG 

WRITE  (2,  920).  1AVG<  J.  1  ).  IAVGI  J,  2).  1AVG1J.  3),  CAVQ1J,  1  ) ,  GAVG  ( J,  7). 
CGAVG<  J.  2) .  GAVC  (  J.  3) .  GAVG  (  J.  4  ) ,  I AVG  (  J.  4  ) ,  GAVG  (  J.  5) .  GAVG  <  J.  6) 

920  FORMAT  (  '  ' .  1 3.  '-  ' .  1 3.  6X,  I  6,  5X.F10.  2.  4X.  F10.  2.  6X.  F7.  2.  6X.  F7.  2. 

C5X.  F7.  2.  SX.  110,  3X.  F12.  2.  4X.  F7.  2) 

930  CONTINUE 

WRITE(2.  910).  12 

CL0SE(UN1T*=2.  DISPOSE- 'SAVE')  1CL0SE  'AVERAGE'  FILE 
C 

999  CALL  SCROLL ( 1 . 24 ) 

STATUS=1  ! RESET  STATUS  TO  USER-INPUT 

CALL  WRTINF  (WRITE  INFO  TO  SCREEN 

CALL  CLRL1N1 18. 23) 

C 

RETURN 

END 
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c 


c 

c. .  . 

c 


451 

c 

c. .  . 
c. .  . 
c. .  . 
c. .  . 

c 


90 


100 


110 

140 

C 

C.  .  . 
C 


C 

C.  .  . 
C.  .  . 
c 


c. .  . 


100 

19Q 

C 

200 


C.  .  . 


210 

220 


SUBROUT  1 NE  CONTIC 
INCLUDE  'EDC_MSARCOM  FDR' 

BYTE  OUTF I L  <  22  ) .  L0W<3),  UPR<3).  EXT  (3) 

INTECER»2  START (25b) ,  STOP!  25b  >.  1  AD JPT ( 5 12 > ,  DROP ( 65536) ,  MAP  (6553b! 
DATA  £XT/'M'.  'C'.  'S'/ 

UN10UE-1 

1DR0P=0 

LOOP  THROUOH  ALL  LINES 

DO  450  1 ROW  = 1 RYTOP .  1RYB0T 
TYPE  451.  1  ROW 
FORMAT ( '  ROW  -  '17) 

SEARCH  LINE  FOR  START  AND  STOP  OF  STRUCTURES  ON  LINE.  THE  ARRAYS 
START(*>  AND  STOP!*)  CONTAIN  THE  STARTING  AND  STOPPING  LOCATIONS 
OF  ALL  the  STRUCTURES  ON  THAT  LINE.  'ISNUrt'  CONTAINS  THE  NUMBER  OF 
STRUCTURES. 


ISNUM=0 

JCOL=IRXTOP 

I F ! THRESH ( JCOL.  IRON)  .  EQ. 
ISNUM=ISNUM+1 
START! ISNUMJ-JCOL 
UC0L=JC0L+1 

1F1JC0L  .  GT.  1RXB0T!  GOTO 
I F (  THRESH (  JCOL.  IRON)  .  EQ. 
STOP ( I SNUM ) *JCQL-1 
JCOL-JCOL+i 

IF (JCOL  .  LE.  IRXBOT)  GOTO 


0)  GOTO  140 


no 

1)  GOTO  100 


90 


SEARCH  EACH  STRUCTURE  FOR  CONTIGUOUS  STRUCTURES  ABOVE 
DO  400  ISTRCT  = 1 .  1 SNUM 

IF1IR0W  .  EQ.  IRYTOP )  GOTO  360  ! TOP  ROW 


LOOP  THROUGH  STRUCTURE.  SEARCH  FOR  ADJACENT  STRUCTURES.  IF  FOUND. 
MAKE  NOTE  OF  UNIQUE  NUMBER  OF  ADJACENT  STRUCTURE. 


I ADJ=0 

DO  300  JCDL=ST ART (ISTRCT). STOP (ISTRCT) 

IF (  THRE5H(  JCOL.  1R0W-1)  .  EQ.  0)  COTO  200 
ADJACENT  POINT  FOUND 
IF  (  I ADJ  .  EQ.  0)  COTO  190 
DO  ISO  J»l.  IADJ 

IF  ( THRESH ( JCOL.  IROW-1)  .  EQ.  IADJPT(  J)  )  OOTO  200  .'UNIQUE  *  ALREADY  FOUND 

CONTINUE 

IA0J=1A0J+1 

I ADJPT ( I ADJ) “THRESH (JCOL. IROW-1 1 


IF ( I ADJCT  .  EQ.  0)  GOTO  300  • DONT  CHECK  FOR  DIAGONALS 

I F ( JCOL-1  . LT.  IRXTOP )  COTO  230  ! DDNT  CHECK  OUTSIDE  ROI 

IF(THRESH( JCOL-1.  IROW-1)  .  EQ.  0>  COTO  230 

ADJACENT  POINT  FOUND 

IF< IADJ  .  EQ.  0)  COTO  220 

DO  210  J-I.  IADJ 

IF(THRESH( JCOL-1.  IROW-1 )  .  EQ.  IADJPT(J))  OOTO  230  IUNIQUE  *  ALREADY  FOUND 

CONTINUE 

IADJ* IADJ* 1 
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c 

230 
C.  .  . 


240 

250 

300 

C 


310 

C 


315 

C 


320 

C 


330 

34  0 
350 

C 

C.  .  . 
C 

360 

370 


C 

C.  .  . 


1 ADJPT! 1ADJ) -THRESH  <  JC0L-1.  1RDW-1  ) 

IF ( JC0L  + 1  .  CT.  1RXD0T)  GOTO  300  ' DONT  CHECK  OUTSIDE  R01 

IF (THRESH! JCOL+1 »  IRQW-1)  .  EQ.  01  GOTO  300 

ADJACENT  POINT  FOUND 

IF (  1 ADJ  .  EQ.  0)  GOTO  250 

DO  240  J=l.  IADJ 

IF(THRESH(JCOL-fl.  1R0W-1)  .  EQ.  2ADJPT(J))  GOTO  300  IUNIQUE  M  ALREADY  FOUND 

CONTINUE 
I  ADJ=1  ADJ-M 

1ADJPT( 1 ADJ) -THRESH  1 JC0L*1.  1R0W-1 ) 

CONTINUE 

I F ( 1 ADJ  . EQ  0)  GOTO  360  'NO  ADJACENT  POINTS.  ASSIGN  NEXT  VALUE 

IF! IADJ  .  EQ.  1)  GOTO  3B0  !  ONE  ADJACENT  POINT 

MORE  THAN  ONE  STRUCTURE  IS  ADJACENT.  SET  PRESENT  STRUCTURE  TO  LOWEST 
ADJACENT  STRUCTURE  VALUE,  THEN  SEARCH  FOR  AND  RESET  ALL  OTHER  PREVIOUS 
STRUCTURE  VALUES  TO  THAT  SAME  LOW  VALUE.  SEARCH  DY  GOING  FROM  PREVIOUS 
LINE  TO  THE  FIRST  LINE  (DOTTOli  TO  TOP).  WHEN  STRUCTURE  VALUE  IS  NO 
LONCER  ENCOUNTERED,  DROP  IT  FROM  THE  SEARCH 


I M I N=999?9 

DO  310  J=l, IADJ  ! SELECT  MINIMUM  STRUCTURE  VALUE 

IF ( I ADJPT ( J )  .  LT.  IM1N)  1M1N-1ADJPT! J) 

CONTINUE 


DO  315  J=l. IADJ  (MAKE  NOTE  OF  DROPPED  VALUES 

IF < I ADJPT ( J)  .  EQ.  1M1N)  GOTO  315 

I  DR OP 3 1  DROP* 1 

DRDP ( 1DR0P > =1 ADJPT (J) 

CONTINUE 


DO  320  JCOL*=ST  ART  (  ISTRCT )  •  STOP  (  ISTRCT  )  '.ASSIGN  MIN  VALUE  TO  'THRESH' 

THRESH ( JCOL. IR0W)«IM1N 

CONTINUE 


DO  350  K=l. IADJ 

IF ( 3  ADJPT  <  K )  .  EQ.  1MIN)  COTO  340  ! DON 'T  DROP  IMIN  VALUE 

DO  340  I«1R0W-1,  IRYTOP.  -1 

IFLC-0 

DO  330  J=IRXTOP,  IRXBOT 

IF ( THRESH! J.  I )  .  NE.  IADJPTIK))  COTO  330 

THRESH! J. I >»IM1N  ! RESET  VALUE 

IFLG=1 

CONTINUE 

IF! IFLG  .EQ.  0)  COTO  350  !  DIDN'T  FIND  ON  LINE.  STOP .SEARCH 

CONTINUE 

CONTINUE 

COTO  400 


ASSIGN  THE  STRUCTURE  NEXT  UNIQUE  VALUES.  NO  ADJACENT  POINTS 


DO  370  JCOL=START( ISTRCT).  STOP(ISTRCT) 

THRESH ( JCOL, IR0W1-UN3QUE 

CONTINUE 

UNIQUE=UNIQUE+1 

COTO  400 


ASS  I CN  THE  STRUCTURE  THE  VALUE  OF  THE  SINCLE  ADJACENT  STRUCTURE 
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c 

380 

390 

C 

400 

450 

C 

C.  . 
C.  . 

c. . 
c.  . 
c.  . 
c 

c 

452 

454 

C 


800 

C 

810 

C 


820 

C 


840 

C 

C.  .  . 
C 


850 
C.  .  . 

870 
C.  .  . 


DO  390  JCOL=START I 1STRCT >. STOP ( ISTRCT) 

THRESH!  JCOL.  IRQW)=1ADJPT( 1 ) 

CONTINUE 
GOTO  400 

CONTINUE 

CONTINUE 

AT  THIS  POINT,  THRESH  CONTAINS  UNIQUELY  LADELLED  MICROSTRUCTURES. 
HOWEVER.  THE  UNIQUE  NUMBERS  MUST  BE  COMPRESSED  SINCE  MANY  VALUES 
WERE  DROPPED.  THE  VALUES  DROPPED  ARE  IN  THE  ARRAY  ' DROr ( I  DROP ) ' . 
THIS  IS  ACCOMPLISHED  BY  CREATED  A  'MAP',  COMPRESSING  THE  MAP  AND 
THEN  USING  THAT  COMPRESSED  MAP  TO  ALTER  THE  VALUES  IN  'THRESH' 

UNKJUE  =  UNIQUE-1 

TYPE  452,  UNIQUE 

FORMAT (  '  NUMBER  OP  MST*  BEFORE  COMPRESSION  -  '110//'*:  'I 

type  454.  unique,  idrop 

FORMAT (  '  UNIQUE  -  '17'  IDROP  -  '  17  > 

ICTMAP^O 

MAP (0)e0  ! CLEAR  MAP 

DO  800  U»l. UNIQUE 
MAP ( J ) -0 

DO  810  UD=1. IDROP  ISETUP  DROPPED  POINTS 

MAP  (  DROP  (UD)  >  —  1 

DO  820  J«l, UNIQUE  !FILL  IN  MAP 

IF ( MAP  I J )  .  EQ  -1)  THEN 
MAP  (  U  )  *=0 
GOTO  820 

END  IF 

ICTMAP=1CTMAP+1 

MAP(J)=ICTMAP 

CONTINUE 

UNIQUE-ICTMAP 

DO  840  IROW=IRYTOP. IRYBOT  IAPPLY  MAP  TO  'THRESH' 

DO  640  JCOL=IRXTOP.  IRXBOT 
1F(THRESH(UC0L.  IROW)  .  EQ.  0)  COTD  640 
THRESH!  JCOL.  1R0W)*>MAP  (  THRESH  (  UCOL.  IROW)) 

CONTINUE 

WRITE  TO  OUTPUT  FILE 

CALL  CONFNM ( OUTF 1L, EXT)  ! GET  OUTPUT  FILENAME 

0PEN(UNIT»3,  NAME-OUTFIL.  TYPE* 'NEW ' ) 

NCHANS* 1  ! ONE  CHANNEL  OF  DATA 

NHEADR-2  ! NUMBER  OF  HEADER  LINES 

WRITE  (3.  650),  NROWS,  NCOLS.  NCHANS,  NHEADR 
FORMAT (414) 

TYPE  CHANNEL  USED.  ROI  INFO,  ADJACENCY  FLAC  AND  THRESHOLD  RANGE 
WRITEO.  670).  ICHAN,  1RXT0P.  IRYTOP.  IRXBOT,  IRYBOT.  IADJCT,  ITLOW.  ITUPR 
FORMAT (BI 7) 

WRITE  'THRESH'  ARRAY  TO  FILE 
DO  603  1-IRYTOP,  IRYBOT 

WRITEO.  680 ) .  (THRESH!  J,  1  ),  J-IRXTOP,  IRXBOT) 


n  o  o  o  o  o  o 


690 

685 
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FORMAT (2216) 

CONTINUE 

CLOSE (UN IT *3.  DISPOSE® 'SAVE ' ) 


RETURN 

END 


o  n  o  o  o  o 
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SUBROUTINE  STAT I S  <  NUMAVG ) 

INCLUDE  'EDC_MSARCOM.  FOR  ' 

INTEGER *2  COORD <2.  8) ,  PERMCT.  ORTHO,  DIAC,  ODma i ,  PERMC 1 .  ICOUNT(S) 

INTEGER  ODma i S 

BYTE  DUTF IL ( 22 ) »  EXT !3 ) 

DATA  COORD/- 1.  -1,  -1.0.  -1.  1.0.  X.  1.  1.  1.0,  1.  -1.0.  -1/ 

DATA  EXT/  'M'.  'S'.  'T'/ 

SQRTW0=2*«0.  5 
C 

C.  .  .  OPEN  FILE  FOR  Ml CRDSTRUCTURE  STATISTICS.  WRITE  HEADER  INFO 
C 

CALL  CONFNM ( OUTFIL. EXT)  'CET  OUTPUT  FILENAME 

OPEN ( UNI T=3. NAME=0UTF1L, TYPE-'NEW') 

WRITEO,  10).  OUTFIL,  ITLOW,  ITUPR 

10  FORMAT ( '  FILENAME  =  '22A1.T40'  LOWER  THRESHOLD  -  '13'  UPPER 

C  THRESHOLD  =  '13) 

WRITE  I  3.  12),  UNIQUE  'HMST.PIX  SIZE.  SIZE  FLAG 


12  FORMAT ( '  '15'  UNIQUE  MICROSTRUCTURES  FOUND.') 

IF! IPIXFG  .  NE.  1  )  THEN 


14 

WRITE!  3.  14).SIZP1X,  IPIXFG 
FORMAT!'  ACTUAL  PIXEL  SIZE  - 

'F7.  3' 

MICRONS. 

FLAG  -  ' 12 ) 

16 

ELSE 

WRITEO,  16).  SIZPIX,  IPIXFG 
FORMAT!'  DEFAULT  PIXEL  SIZE  - 

‘F7.  3  ' 

MICRONS. 

FLAG  -  '12) 

22 

END  IF 

WRITEO.  22) 

WRITE! 3,  24) 

WRITEO.  26) 

FORMAT!/// ' 

WEIGHTED 

WEICHTED 

C  MAXIMUM  MAXIMUM') 

24  FORMAT!'  UNIQUE  X  Y  X  Y 

C  PERIMETER  PERIMETER  HORIZONTAL  VERTICAL 


C  AVERAGE') 

26  FORMAT!'  MICHOSTRUCTURE  CENTER  CENTER  CENTER  CENTER 

C  AREA  COUNT  LENGTH  RUN-LENGTH  RUN-LENGTH  IOD 

C  IOD  ODma  x ' / ) 

C 

C.  .  .  CLEAR  SOME  VALUES 
C 

IAREAS=0 
XCEN1=0 
YCEN1-0 
XWCEN1-0 
YWCEN 1  *=0 
XTCENT=0 
YTCENT-0 
PERMC 1 *0 
PERIM1-0. 

IMAX_X_SUM=0 
IMAX~Y_SUM=0 
10DSUH«=0 
ODmaiS-O 

CALCULATE  STATISTICS  FOR  ALL  MICROSTRUCTURES 
DO  900  IMST«1.  UNIQUE 

SEARCH  FOR  FIRST  OCCURANCE  OF  MICROSTRUCTURE  #IMST 
DO  100  1R0W=IRYT0P.  IRYBOT 


n  o  n  o  o 
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100 

c 

c. . 
c. . 
c. . 
c. . 
c. . 
c. . 
c 

c. . 
c. . 
c. . 
c. . 
c 
c 

1 10 


c 

120 


125 


c. .  . 


c 

160 


DO  100  1C0L=IRXT0P.  1RXB0T 

I F ( THRESH ( I  COL.  IRQW)  .  E<3.  IMST)  GOTO  110 

CONTINUE 

CALCULATE  PERIMETER  BE  SEARCHING  ADJACENT  POINTS.  IF  1ADJCT-0. 
SEARCH  112,4,6,8.  IF  1ADJCT-1.  SEARCH  HI.  2,  3,  4,  5.  6.  7,  8. 

1  2  3 
B  •  4 
7  6  5 

SINCE  FIRST  POINT  WAS  FOUND  BY  SEARCHING  FROM  THE  LEFT  (FROM  #B) 
UITH  NOTHING  IN  THE  PREVIOUS  LINE.  THE  FIRST  SEARCH  POINT  IS 
K4.  AFTER  A  POINT  IS  FOUND,  SEARCH  FOR  THE  NEXT  POINT  BY  STARTING 
AT  (PREVIOUS  SEARCH  #  ♦  4  +1 )  OR  <P.  S.  »  -  4  ♦  1  )  IF  P.  S.  d  >■  4. 


IPSRC=4 
I SRC= I PSRC 
PERMCT=0 
0RTH0=0 
DIAC=0 
1 YM I N= I  ROW 
1 YMAX=1R0U 

ixm;n=icol 

I XMA  X  —  I  COL 
IYST=1R0U 
IXST-1C0L 

I F  (  I YST+COORD  (  1 ,  ISRC)  .  LT.  1RYT0P  .OR.  I  YST +CDDRD  (  1 , 

CISRC )  .  GT.  IRYBOT)  GOTO  125  (CHECK  IF  DUT  OF  ROI 

IFUXST+C00RD<2.  ISRC)  .  LT.  IRXTOP  .OR.  I  XST+CDORD  ( 2, 

CISRC)  .  GT.  I RXBOT )  GOTO  125 

IF  (THRESH l  1  XST «-COQRD  12,  ISRC  ) ,  1V6T+CQDRD ( 1 ,  ISRC  >  )  .  NE.  0)  GOTO  160 
ISRC=ISRC+1 

I F ( I  AD JCT  .  EQ.  0)  ISRC-1SRC+1  ! ONLY  CHECK  ORTHO  POINTS 

1F(ISRC  . GT.  B)  ISRC-ISRC-B 

IF ( I  SRC  .  NE.  I PSRC )  GOTO  120 

POINT  IS  ISOLATED.  PERIMETER  IS  1  POINT 

PERMCT- 1 

GOTO  200 

IF(JMOD( ISRC.  2)  .  EO.  0)  THEN  ! DETERMINE  IF  ORTHO  OR  DIAO 

ORTHQ=ORTHO+t 

ELSE 

DIAC-DIAG*! 

END  IF 

IYST  =  IYST*COORD( 1,  ISRC)  IMOVE  TO  NEW  POINT 

I XST- I XST+COORD ( 2, ISRC) 

IF  ( THRESH(  J  XST.  IYST)  .  GT.  0)  THEN 

PERMCT=PERMCT*1  ! COUNT  'NEW'  POINT 

THRESH! 1XST,  IYST)-  —THRESH! I XST.  IYST)  INECATE  POINT  AS  FLAO 

END  IF 

AT  THIS  POINT,  YOU  MUST  CHECK  IF  ALL  SURROUNDING  POINTS  HAVE  BEEN 
COUNTED  SEARCH  AND  SEE  IF  ALL  SURROUNDING  POINT  OF  THE  PRESENT 
MICROSTRUCTURE  HAVE  BEEN  SET  TO  -VALUE 


! COUNT  OF  PERIMETER 
(COUNT  OF  ORTHOGONAL  POINTS 
(COUNT  OF  DIAGONAL  POINTS 

(set  parameters  for  blob  range 


IFUYST  .  EQ.  IRON  .AND.  IXST  .  EQ.  ICOL)  THEN 


ooonoononnnonononno 
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1 ADJ_C0UNT=0 

IADU_N0T_C0UNT«0 

ISKIP=2 

IFUADUCT  ,  EQ.  1)  ISKIP=1 
DO  170  ICHK=ISK1P. 8, 1SKIP 

IF! I YST*COORD<  1.  ICHK)  .  LT.  IRYTOP  .OR.  I YST+COORD  ( 1 , 

CICHK)  . CT.  1RYB0T)  GOTO  170  ! CHECK  IF  OUT  OF  ROI 

IF<  IXST+C00RD<2.  ICHK)  .  LT.  IRXTOP  .OR.  IXST+C00RDI2, 

CICHK)  .  GT.  IRXBOT)  GOTO  170 

1FITHRESH! IXST+C00RDI2.  ICHK).  IYST+COQRD! 1.  ICHK) )  .  LT.  0)  THEN 
I  ADU_COUNT  =  I  AD  J_C0UNT  +  1 
ICOUNT! 1ADU_C0UNT)-ICHK 

END  IF 

IF(THRESH( IXST+COORDI2.  ICHK).  I YST  +  COORD ! 1  ■  ICHK) )  .  CT.  0) 

C  I ADU_NOT_COUNT=I ADJ_NQT_C0UNT+1 
170  CONTINUE 

IF( I ADJ_NOT_COUNT  .  EQ.  0)  GOTO  200  (BACK  AT  1ST  POINT 

IF! 1 ADJ~COUNT  .  EQ.  1)  GOTO  180  ! CONTINUE  SEARCH 

DO  172  UU-1 .  IADJ_C0UNT-1 
DO  172  KK*2. IADJ_COUNT 
1D1F«=IADS<  ICOUNT  l  JJ )  -  I  COUNT  i  KK )  ) 

IF(ID1F  .  EQ.  1  .OR.  IDIF  .  EQ.  7)  GOTO  172 

GOTO  200  ! SEPERATED  BY  TWO.  BACK  AT  1ST  POINT 

172  CONTINUE 

GOTO  180 
END  IF 


ISK1P  =  2 

I F  <  1  ADUCT  .  EQ  1  )  1 SK I  P*=  1 
DO  170  ICHK*=ISKIP.  8.  ISK1P 
I F ( I YST  .  LT.  60)  GOTO  173 

TYPE  172.  1 XST .  IYST.  ICHK,  THRESH! IXST+C00RDI2,  ICHK).  I YST+COORD! 1 .  ICHK)  > 

172  FORMAT!'  I XST  -  '15'  IYST  -  '15'  ICHK  -  '14'  THRESH  -  'I7> 

173  IFUHRESH!  IXST+C00RD12.  ICHK).  IY6T+C0DRD! 1  ■  ICHK))  .  EQ.  0)  THEN 

IZRFLC«0 
GOT □  170 

END  IF 

IF ! THRESH! I XST +COORD ! 2.  I CHK) ,  I YST+COORD ! 1 ,  ICHK >  >  .  LT.  0)  IZRFLC-t 
IFITHRESH! IXST+C00RDI2.  ICHK),  IYST+CQORDI 1.  ICHK))  .  GT.  0 
C.  AND.  IZRFLG  .  EQ.  0)  COTO  180 
170  CONTINUE 

GOTO  200  '.BACK  AT  1ST  POINT 

END  IF 


C. . .  SET  RANGE  VALUES 

IBO  IF ( IYST  .LT.  IYMIN)  1YM1N-1Y5T 

IF! IYST  .  GT.  IYMAX)  IYMAX-IYST 
IF! I XST  .  LT.  IXMIN)  IXMJN»IXST 
IF < I XST  . GT.  IXMAX)  1XMAX-IXST 
C.  .  .  CALCULATE  WHERE  TO  LOOK  NEXT 
IF! 1SRC  .  LE.  A)  THEN 
IPSRC-ISRC+4 

ELSE 

IPSRC-ISRC-A 

END  IF 

IPSRC=IPSRC+2 

1FUPSRC  .  CT.  B>  1PSRC-1PSRC-B 
ISRC-IPSRC 
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c 

c. . 
c 

200 

c 

c. . 
c. . . 

c 


290 


295 

300 


320 

C 

C.  .  . 
C.  .  . 
C 

C.  .  . 


C 


390 


C0T0  120 

CALCULATE  PERIMETER  DISTANCE 

PER1M«(ORTHO+<DIAC»SQRTWO>  )*SIZP1:'  1N0RMALIZED  TO  MICRONS 

SECOND  PASS  THROUGH  RANGE  THAT  CONTAINS  MICROSTRUCTURE. 

DETERMINE  MAX  HORIZONTAL  AND  VERIT1CAL  DIAMETERS 

IMAX_Y=0 

IMAX~X*=0 

DO  320  1C0L-IXMIN,  IXMAX  'CALCULATE  Y-MAX  DIAMETER 

IST=0 

DO  300  IROW=IYMIN,  IYMAX 

IF( IADS(THRESH( ICOL.  IRON) )  .  EQ.  IMST)  THEN 
1F<  1ST  .  EQ.  1  )  GOTO  290 
IST=1 

I STLC« I ROW 
CONTINUE 

ELSE 

I F  t  1ST  .  EQ.  0)  GOTO  295 

IF< ( IROW-1 >-3STLC  .  CT.  IMAX_Y>  IMAX  Y-1R0U-1-1STLC 

IST=*0 

CONTINUE 

END  IF 
CONTINUE 

IF<  1ST  .  EQ.  0)  GOTO  320 

IF( ( IROW-1 I-1STLC  .GT.  1MAX_Y>  IMAX  Y-IROU-l-ISTLC 

CONTINUE 

THIRD  PASS  THROUGH  RANGE  THAT  CONTAINS  MICROSTRUCTURE.  DETERMINE 
GEOMETRIC  CENTER,  AREA,  1UD.  ODman,  AND  AVERACE  OPTICAL  DENSITY. 

CLEAR  ALL  ARRAYS 

IAREA=0 

IYSUM=0 

IXSUM=0 

WYSUM=0. 

UXSUM-O. 

IOD-O 
ODma i«-l 


DO  420  IROW-IYMIN,  IYMAX 
IST-0 

DO  400  IC0L-1XM1N,  IXMAX 

IFUABSITHRESHUCOL,  IROWM  .  EQ.  IMST>  THEN 
I  IMG” IMAGE ( I  COL,  IROW) 

IF  <  1ST  .  EQ.  1 )  GOTO  390 
IST-1 

ISTLC-ICOL 
IAREA-IAREA+1 
IYSUM«IYSUM»IROW 
IXSUM-IXSUM*ICOL 
IOD=IOD*I IMG 

HYSUM»WYSUM+FLOAT< IROW*IIMC> 

WXSUM»UXSUMfFLOAT< ICDL*IIMC> 

I F ( I  IMG  . CT.  ODmai)  ODmai-IIMC 

ELSE 

1FUST  .  EQ.  0)  GOTO  395 

I F  <  (  ICDL-1 )— ISTLC  .CT.  IMAX_X)  IMAX  X-ICOL-1- .  STLC 
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395 

400 

420 


C 

C.  .  . 
C 


eoo 

c 

c. . . 
c. . . 
c 


020 

030 


C 

900 

C 

C.  .  . 

c 


IST  =  0 
CONTINUE 

END  IF 
CONTINUE 

IF< 1ST  .  EQ.  0)  GOTO  420 

IF<  ( ICOL-1  I-ISTLC  ,CT.  1MAX_X) 

CONTINUE 

IMAX_Y«IMAX_Y*1 

IMAX~X=IMAX_X+1 

YCENT=FLOAT( 1 YSUM ) /FLOAT < I AREA ) 

XCENT =FLOAT ( 1 X SUM) /FLOAT  f  2  AREA) 
YWCENT=WYSUM/FLOAT( IOD) 
XUCENT-WXSUM/FLOAT i IOD) 

AREA-FLOAT  (1AREA)«(£I2PIX«SI2PIX> 
AVCOD-FLOAT  < IOD) /FLOAT ( 1  AREA) 

IF ( 1  AREA  .  EQ.  PEKMCT)  THEN 
PERMCT--1 
PER1M— 1.  0 

END  IF 


IMAX  X-I COL— 1 - 1 STLC 


!  GEOMETRIC  CENTERS 


.'NORMALIZED  AREA 
( AVERAGE  OPTICAL  DENSITY 
(NO  'INTERIOR'  POINTS 


WRITE  INFORMATION  TO  FILE 


WRITEI3.  BOO),  IMST.  XCENT.  YCENT,  XUCENT.  YWCENT,  AREA,  PERMCT, 
CPERIM,  IMAX_X,  IMAX_Y.  IOD, AVGOD.  ODmai 

FORMAT  (  '  ' T6,  15,  T10.  F7.  2.  T27,  F 7.  2,  T36,  F7.  2,  T45.  F7.  2,  T53,  F9.  2, 
CT64.  15,  T74.  F0.  2,  T89.  13.  T102,  13,  TI09,  10,  T118,  F7.  2,  Tl27,  23) 


KEEP  A  RUNNING  TOTAL  09  ALL  STATISTICS  FDR  THIS  THRESHOLD  RANGE. 

THIS  WILL  BE  WRITTEN  TO  AN  OUTPUT  FILE  AT  COMPLETION  OF  ALL  PARTICLES 


XCEN1 -XCEN1 ♦  XCENT 

YC£N1=YCENJ+YCENT 

XTCENT  =  XTCENT-MXCENT»10D> 

YTCENT-YTCENT+< YCENT* IOD) 

XWCEN1-XUCEN1+XWCENT 

YUCEN1 =YWCENI ♦YWCENT 

1  ARE AS* I  ARE AS ♦ I AREA 

IF  (PERMCT  .  LT.  0)  GOTO  020 

PERMC I -PERMC I +PERMCT 

IFIPERIM  .  LT.  0)  GOTO  030 

PERIM1-PER1M1+PERIM 

IMAX_X_SUM»IMAX_X_SUM+IMAX_X 

IMAX~Y_SUM=IMAX~Y~SUM+lMAX~y 

iodsum»iodsum+iod“ 

QDma  «S«ODma  « S+ODma < 

CONTINUE 

CLOSE (UNIT-3. DISPOSE-  'SAVE  '  > 


•X  AND  Y  CENTERS 

(WEIGHTED  X  AND  Y  CENTERS 
(TOTAL  AREA 

(TOTAL  PERIMETER  COUNT 

(TOTAL  PERIMETER  LENGTH 
(TOTAL  X,  Y— RUN  LENGTH 

(TOTAL  IOD 
(TOTAL  ODmai 


CALCULATE  'AVERAGE'  OF  STATISTICS  AND  WRITE  70  UNIT  2 


IF( UNIQUE  .  LT.  1)  THEN 
AUNO-O 
AOUTl-O 
A0UT2-0 
A0UT3-D 
A0UT6-0 
A0UT7-0 
GOTO  950 


END  IF 


o  o  o  n 
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AUNQ'FLOAT  <  UN  I OUE ) 

AOUTl'FLOATI 1AREASI/AUNQ 
A0UT2=XCEN1/AUNQ 
AOUT3=YCENl /AUNQ 
AOUTt'XTCENT/FLOATI  10DSUM) 

AOUT7=YTCENT/FLOAT ( IODSUM) 

950  WRITE  (2.  920) ,  1TL0W.  ITUPR.  I  AREAS,  UNIQUE.  ADUT1,  A0UT2.  A0UT3, 

CA0UT6,  A0UT7 

920  FORMAT  I  '  '.  13.  13.  6X,  16.  11X.  16.  11X.  F9.  2,  5X,  F7.  2.  6X,  F7.  2. 

C6X,  F7.  2.  7X,  F7.  2) 

C.  .  .  STORE  STATS  FDR  SECOND  SECTION  OF  OUTPUT 

NUMAVC'NUMAVG-*- 1 
lAVCINUMAVC,  1  )  =  1TL0U 
I AVG ( NUMAUC.  2 ) “ ITUPR 
IF  t UNI OUE  .  LT.  1)  THEN 

IAUGINUMAVG.  3)=0 
I  A  VG  I NUMAVG .  41=0 
GAUG  (  NUMAUG,  1  )'0. 

G AUG I NUMAVG .  2>=0. 

GAVGINUMAVG,  3>'0. 

GAVGINUMAVG.  4  >'0. 

GAUG  (  NUMAVG.  5  )  *=0. 

G AUG < NUMAVG , 61=0. 

GAUG  (NUMAVG,  7  >'0. 

ELSE 

IAUGINUMAVG.  3 ) =PERMC 1 

IAUGINUMAVG.  4)  =  I0DSUM 

GAUGINUMAUG.  1 > 'FLOAT  I PERMC 1 ) /AUNQ 

GAUG < NUMAVG,  2)=PERIM1 /AUNQ 

GAVG I NUMAUG.  3) “FLOaT ( 1 MAX_X_SUM ) /AUNQ 

GAUG I NUMAUG , 4 ) 'FLOAT I 1MAX_Y_SUM ) /AUNQ 

GAVG I NUMAUG.  5 ) 'FLOAT  I  I ODSUM ) /AUNQ 

GAUG I NUMAVG. 6 ) 'FLOAT  I ODmaiS) /AUNQ 

GAVCINUMAVG. 7)"PER1M1 

END  IF 
999  RETURN 

END 
C 


uuuuu 
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c 


700 

710 

3 


715 


C 


SUBROUTINE  CONFNM ( OUTF I L. EXT ) 

INCLUDE  'EDC_MSARCOM.  FOR  ' 

BYTE  0UTF1L122).  LOW l 3 > .  UPR<3><  EXT13) 


DO  700  JP01NT=1 .  14 
1F(FNAME<  JP01NT)  .  EQ.  ') 
OUTF IL ( UP01NT )“FNAME<UPOINT 


COT 0  710 
) 


CONTINUE 

ENCODE  3. 3. LOW) .  ITLOW 
ENCODE  (  3,  3/  UPR  ) .  1  TUPR 
FORMAT  t 13) 

0UTF1L<UP0INT>= 

OUTFIL( JP Q I NT*4)~ 

DO  715  J=l.  3 

1 F  (  LOW  (  J  )  .  EQ.  '  ')  L0W(J)  =  '0' 

IF(UPR<U)  EQ.  ‘  ')  UPR(J)“'0' 

OUTF 1 L  <  JPOINT* J ) =LDW ! J ) 

OUTF IL ( JPOJNT+J+4 ) =UPR ( J) 


CONTINUE 

OUTFIL(UPO1NT*0)=  ' 
OUTFIL<UPOINT*9>-£XT< ’ ) 

OUTF IL 1  UP 01 NT* 1 0 ) =EXT ( 2) 
0UTFIL(UP01NT*1  1  )<=EXT  (3) 


'COPY  FILENAME 


return 

END 
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SUBROUTINE  WRT1NF 
INCLUDE  'EDC_MSARCOM,  FOR  ' 

BYTE  BOLD*  4) .  ATTOFF  *  4 ) , REVERSC4).  BLINK <4 ) 

DATA  BOLD/27,  91. 49.  109/  (ESC  t  I  m 

DATA  ATTGFF/27. 91, 48. 109/  ! ESC  C  0  m 

DATA  REVERS/27,  91.  55,  109/  (ESC  C  7  m 

DATA  BLINK/27.91,53,109/  !ESC  C  5  m 

IT0P=4 

CALL  VTCURP < ITOP, 1 1 ) 

TYPE  100.FNAME 
CALL  VTCURP  (  I  T0P->-4,  10  > 

TYPE  140.  ITLOW 
CALL  VTCURP* IT OP +  4,  27) 

TYPE  140. 1TUPR 
CALL  VTCURP ( I TOP+6.  15) 

TYPE  140. ITLINC 
CALL  VTCURP ( I TOP  +  6,  32 ) 

TYPE  140. ITUINC 
CALL  VTCURP * ITOP  +  4,  48) 

TYPE  140. IRROW 
CALL  VTCURP* ITOP  +  4,  69) 

TYPE  140.  1RC0L 
CALL  VTCURP* IT 0P  +  5.  49 ) 

TYPE  140.  I R  XT  OP 
CALL  VTCURP* ITOP  +  S, 67) 

TYPE  140.  1RYT0P 
CALL  VTCURP* ITOP+6.  52) 

TYPE  140, IRXBOT 
CALL  VTCURP* ITOP  +  6,  70) 

TYPE  140, IRYBOT 

CALL  VTCURP* IT OP +  5.  22)  ! THRESHOLD  INCREMENT 

IF  ( I  T 1  NCR  .  EQ.  1  )  THEN 
TYPE  150 

ELSE 

TYPE  160 

END  IF 

■'LL  VTCURP* ITOP  +  8, 21 )  ! ADJACENCY  CRITERIA 

jF< IADJCT  EO  1 )  THEN 
TYPE  170 

ELSE 

TYFE  ISO 

END  IF 

CALL  VTCURP* ITOP+9,  18) 

IF* IP1XFC  .  EQ.  1 )  THEN 
TYPE  i 90 

ELSE 

TYPE  200, SI2PIX 

END  IF 

IF ( STATUS  .  EQ.  1)  THEN  ‘PROCRAM  STATUS 

TYPE  230. ATTOFF 
CALL  VTCURP*  1TQP  +  H.  171 
TYPE  210 


IPIXEL  DIMENSION 

1USER  DID  NOT  KNOW  DIMENSION 


‘FILENAME 
1  LOWER  THRESHOLD 
! UPPER  THRESHOLD 
! LDWER_I NCR 
!UPPER_INCR 
!  ROMS 
!  COLUMNS 
(X-TQP 
! Y-TOP 
! X-BDTTOM 
! Y-BOTTDM 


ELSE 


TYPE  230,  BOLD 
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TYPE  230.  REV'ERS 

TYPE  230. SlINK 

CALL  VTCURP!lTOP+ll.  17> 

TYPE  220 

TYPE  230.  ATTQFF 

END  IF 
C 

100  FORMAT  (  1  4  A 1  ) 

120  FORMAT  (  '*  '■  11) 

i  40  FORMAT  (  13) 

150  FORMAT!  '^-ENABLED  ')  , 

160  FORMAT ! '^DISABLED' ) 

170  FORMAT! '+ORTHOGONAL/D I  AGONAL ' ) 

ISO  FORMAT! ' +0R THDGON AL  ') 

190  FO°MAT ( '+*****♦* ' ) 

200  FORMAT!  F7.  3) 

210  FORMAT ! '♦USER-INPUT ' ) 

220  FORMAT! '^EXECUTING  ') 

230  FORMAT  <  4A1  ) 

C 

C 

RETURN 

END 

C 

C 

c 

c 

c 


noooo  a*-  rv)  non 
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SUBROUTINE  VTCURP ( LINE.  COL ) 

PLACE  CURSOR  AT  LINE  b  COL  OF  SCREEN 
IMPLICIT  INTEGER«2(A-Z> 

BYTE  CP(B)  iHORIZONTAL  b  VERTICAL  CURSOR  POSITION 
DATA  CP/27.  91. 48.  48.  59.  4B.  AS.  72/  'ESC  COO)  0  0  H 
ENCODE  (2.  2.  CP  (3)  ILINE 
ENCODE  <  2,  2.  CP<6>  )COL 
FORMAT (12) 

IF  <CP<3>  .  EO.  *  ')  CP  ( 3 )  «  'O' 

IF  ( CP  <  6 )  .  EQ.  '  ')  CP  (6)  -  'O' 

TYPE  10,  CP 
FORMAT  <  '*  '8A1 ) 

RETURN 

END 
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c 

100 


c 

c 

c 

c 

c 


c 


100 

200 


c 

c 

c 

c 

c 


c 

100 


c 

c 

c 

c 

c 


SUBROUT  1 NE  CLRTXT 
BYTE  COMAND(A) 

DATA  COMAND/27,  91.  50.  74/ 

TYPE  lOO.COMAND 
FORMAT ( '♦'4A1 ) 

RETURN 

END 


SUBROUTINE  CLRL I N { IT  OP ,  1  DOT } 

BYTE  MOVE (  7  > » ERASE ( A ) 

DATA  M0VE/27,91.0,0,59,49,72/ 

DATA  ERASE/27,91,50,75/ 

DO  200  J-ITOP.  IBOT 
1LINE=J 
MOVE ( 3  > “ '0 ' 

IFIILINE  .  CT  9)  M0VE(3)-(  ILINE/101+AB 
MOVE  <  A ) =  JMOD ( 1 L 1 NE,  101+4B 

TYPE  100, MOVE  IMOVE  CURSOR  TO  LINE 

TYPE  100. ERASE  ! ERASE  ALL  OF  LINE 

FORMAT ( '+ '7A1 > 

CONTINUE 

RETURN 

END 


! ESC  C  PI  I  1  H 
! ESC  L  2  K 


SUBROUTINE  CLFSCR 
BYTE  1 NI  T  13),  END  <  2  I 

DATA  INIT/27,  80.  112/  !E5C  P  p 

DATA  END/27,92/  !ESC  \ 

TYPE  100,  INIT.  END 
FORMAT  1  '+'3A1,  'S(E)  2A1) 

RETURN 

END 


n  n  o  o  o 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  VTDBLE <  TEXT.  NCHARS,  LINE,  COL,  VID) 
DOUBLE  WIDTH  VT  TERMINAL  TEXT  PROGRAM 


CALLINC  VARIABLES: 

TEXT  -  input  (tring  of  tut 

NCHARS  -  input  number  of  character*  of  the  teit  C0-A03 

if  nchan  -  O,  (imply  **t  the  eiliting  tot  to  double  width 
LINE  —  input  line  of  teit  C  I  -  23  3 
COL  -  input  (tarting  column  of  teit  t  1  —  AO  3 

VID  —  normal  t0>  or  revert*  il)  video 

IMPLICIT  1NTEGER«2(A-Z ) 

EVTE  DW(3).  TEXT  (401,  CP  ( 8 ) .  R  v  (  A  ),  Nv<  3 ) 

DATA  CP/27,  91.  48,  AB,  59,  AS,  48,  72/  [CURSOR  POSITION 
DATA  DW/27.35.5A/  'DOUBLE  WIDTH  SINGLE  HEIGHT 
DATA  RV/27. 91. 55, 109/  ! ESC  C  7  m  REVERSE  VIDEO 

DATA  NV/27, 91, 109/  ! ESC  t  m  NORMAL  VIDEO 


IF 

(LINE 

.  LT. 

1>  LINE 

-  1  !VT  SCREEN  LINES 

FROM  1 

IF 

(LINE 

.  GT. 

23)  LINE 

-  23  ! 

TO  23 

IF 

(COL 

.  LT. 

1  1  COL 

-  1 

IF 

(COL 

.  CT. 

AO)  COL 

*  AO 

IF 

(NCHARS 

*  COL 

.  GT.  Al) 

NCHARS  -  A1  -  COL 

ENCl)DE<2,  2.  CP<3>  )  LINE 
ENCODE (2, 2. CP ( 6 ) )  COL 
2  FORMAT < 12) 

IF  <CP(3)  .  Ed.  •  CP  1 3 )  -  '0'  ! ZERO  FILL  1  D1CIT  NUMBERS 

IF  JCP  (6)  .  EO.  '  ')  CP  ( 6 )  -  '0' 

IF  (NCHARS  .  LE.  0)  THEN  [mat*  exittlnj  tett  double  width 

TvpE  10,  CP.DW 
.  ‘turn 
end  if 

if  Vi d  .  eq.  1)  then 

type  10.  Cp  ,  Dw,  Rv,  (Te  I  t  (  1 ) ,  1-1,  Nchan  ) ,  Nv  [revere*  video  tot 
•  1 '  e 

type  10, CP, DW, Nv.  <TEXT< I ), 1-1. NCHARS)  [normal  video  teit 

eni  if 

10  FOR  MAT  (  58A1 ) 

RETURN 

ENT 


r>  n  n  n  n  r>  no  o  r>  nnnnnnnono 
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asynchronous  keyboard  input  usinc  system  services  oio 

WITH  WAIT  FOR  EVENT  FLAG  TO  PREVENT  POLLING. 


FUNCTION  IESC  (  CHAR  ) 

IMPLICIT  INT£CER«2  (  A  -  Z  ) 
LOG  I  CAL  *  1  FLAG 
INTEGERS  IESC_EFN 
INTEGERS  IOSdTa) 

1NT£CER«2  1ESC_CHAN 
1 NTEGER «2  NEXTM) 

BYTE  FIRST_CHAR 

BYTE  LAST_TW0_CHARS(2) 

CHARACTERS  TERMINAL  /  'TT'  / 


DATA  FLAG  /  .  FALSE.  / 

INCLUDE  M  »IODEF  )'  ! DEFINE  SYMBOLIC  NAMES  FOR  I/O  FUNCTION 

CODE  1 0*READVBLK  AND  IO*M_NOECHO.  THE 
CODE  MODIFIER. 


10 


C 

C 


1 

2 

0 

4 

5 

6 
7 

e 

? 

l 


CHAR  -  0 
IESC  «  0 
IESC_EFN  =  2 
IF  (FLAG)  GO  TO  10 

1STAT  »  SYSt ASS  I GN  ( TERMINAL.  I ESC_CHAN.  ,  > 

IF  (  .NOT.  ISTAI  I  CALL  LIB*STOP  (  XVAL  (ISTAT  )  ) 

FLAG  *  .  TRUE. 

IO_FUNC  -  IO«_READV0LK  .OR.  10*M_NOECHQ  OR.  10»M_ESCAPE 
10_TEST  -  10*_SENSEM0DE  .OR.  I 0*M_TYPEAHDCNT 
ISTAT  -  SYS*QrOW(XVAL( IESC_EFN)  . “  !  READ  A  CHARACTER 

XVAL( IESC_CHAN).  !  WHEN  KEYED. 

XVAL(IO_FUNC  I.  ! 

xrefiioSb  I.  ! 


XREF ( F IRST_CHAR ) > 
XVAL ( 1 ) . 


.  ) 


1 

2 

3 


IF  (.NOT.  ISTAT  )  CALL  LlBtSTOP  (XVAH  ISTAT))  •  STOP  IF  OIO  NOT  SUCCESS 
CHAR  -  FIRST_CHAR 

ISTAT  -  SYS*QIOW(XVAL( 1ESC_EFN).  ! HOW  MANY  MORE  CHARACTERS 

XVAL< IESC_CHAN). 

XVALI 10_TEST)  , 

XREF (  _IOSB)  . 


r>  n  n  n  n 
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c 

c 


c 


4 

5 

6 


1 

2 

3 

4 

5 
b 


' XREF ( NEXT )<.<<•  ) 

IF(  NEXT(l)  .  EG.  0  )  RETURN 

1STAT  "  SYS*aiOW<>/.UAL<lESC_EFN>. 

Y.UAL<  1ESC_CHAN>, 
XVALI  1ESC_FUNC). 
•/.REF  (  I  OSB  )  / 


! GET  THE  REST 


'  y.REF  ( LAST_TWO_CHARS ) .  ‘/.VAL 1 2  ) .  ,  .  .  ) 


CHAR  **  LAST_TWD_CHARS<2 > 

IESC  -  1 


C  WRITE  (6.100)  IESC. IDUM 
CCIOO  FORMAT  (6X,  'BYTE  COUNT: 
C 

20  RETURN 
END 


.14.  '  FIRST  CHARACTER:  A1  ) 
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SUBROUTINE  SCROLL < I TOP. 1 BOT ) 

SET  SCROLL  REGION  TO  ONE  DEFINED  BY  ITOP, IBOT.  THEN  MOVE 
CURSOR  TO  'HOME'.  TOP  LINE  OF  SCROLL  REGION. 

BYTE  INI T ( 2 )  i  STRING ( 6 ) 

DATA  INIT/27.  91/  !ESC  C 

STRING! 1>«< ITOP/lOJ+48 

STRING (2)-( I  TOP-! (STRING! 1 >-48>*10> >+4B 
STR I NC I  3 ) - '  I  ' 

STRING! 4 >-! I BOT/ 10) +48 

STRING! 5)-! IBOT-I ( STR ING ! 4 ) -4B ) *1 0 ) ) *48 
STRING!*)- 'r ' 

TYPE  10,  "33.  (STRING!  J).  J-i.  b) 

FORMAT!  Al.  ' t  ' .  *A  1  ) 

TYPE  20,  INIT 
FORMAT! '♦ '2A1.  "?6h ' ) 

RETURN 

END 


